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Abstract: We introduce fractional flat space, described by a continuous geometry 
with constant non-integer Hausdorff and spectral dimensions. This is the analogue 
of Euclidean space, but with anomalous scaling and diffusion properties. The basic 
tool is fractional calculus, which is cast in a way convenient for the definition of the 
differential structure, distances, volumes, and symmetries. By an extensive use of 
concepts and techniques of fractal geometry, we clarify the relation between fractional 
calculus and fractals, showing that fractional spaces can be regarded as fractals when 
the ratio of their Hausdorff and spectral dimension is greater than one. All the results 
are analytic and constitute the foundation for field theories living on multi-fractal 
spacetimes, which are presented in a companion paper. 
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|6|. Discussion 



1. Introduction 



1.1 Motivation 



The quest for a quantum theory of gravity has reached such a level of sophistication 
that several independent approaches to this problem have been developed by now. 
Each proposes its own physical picture, but they are often related one to the other 
[[!]]. Apart from string theory, we mention in particular group field theory 0, loop 
quantum gravity (LQG) [§, [|, spin foam models || [|, asymptotically safe gravity 
(or quantum Einstein gravity, QEG) |7[|-[^0[|, and simplicial quantum gravity |TI 



the latter having a particular incarnation in causal dynamical triangulations 
(CDT) HH| . All these approaches, starting from different assumptions and using 
different techniques, have achieved important results concerning both the kinematical 
description of quantum space and its dynamics. However, many formal aspects of 
such dynamics are still to be understood and we still lack a conclusive proof that, 
within any of these models, the proposed quantum dynamics of space as a discrete 
entity lead to a continuum description of the same, and to the general relativistic 
dynamics at large scales and in the classical limit. 

Another important issue is the extraction of effective models of quantum gravity 
dynamics for both spacetime and matter. These effective descriptions should be used 
to predict new phenomena and quantum gravity corrections to known particle and 
astroparticle physics, as well as to large-scale cosmological scenarios. An example of 
the type of phenomenology one might deal with is dimensional flow. It has been no- 
ticed that independent models such as CDT ]T5|, [TJJ], QEG (T7J, and Hofava-Lifshitz 
gravity [T^, all exhibit a running of the spectral dimension of spacetime such 
that at short scales the physics is effectively two-dimensional (see also This 
number is not accidental and plays an important role in the construction and renor- 
malizability properties of quantum gravity pl|-[23|. The change of dimensionality at 



different scales and its acquiring non-integer values is typical of multi-fractals, so it 
is customary to describe dimensional flow as a "fractal property" of spacetime. The 



idea that spacetime is "fractal" in extreme regimes has been hovering for a while |24 
especially in relation with the classical mixmaster behaviour of the BKL big bang 
singularity j2"5|-||32|| and with the notion that renormalization properties of gravity 
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improve on a microscopic "spacetime foam" [33, 34]. The concept of fractal spacetime 
has been often shrouded in a halo of vagueness, devoid of any clear-cut definition, 
and it has begun to be realized concretely in quantum models only in recent times. 
Aside from the above cases, there are hints of fractal spacetime behaviour at high 



energies also in effective K-Minkowski non-commutative field theories These 
findings prompted similar investigations in the context of LQG spin foam models, 
with some preliminary results obtained using certain approximations to the full spin 
foam dynamics |]36|-|38|. 

The scenario is that of a fundamental dynamics where the usual notions of space, 
time and geometry emerge only in specific regimes and approximations of the theory. 
At high energies/small scales, the effective dimension of spacetime is two, while at 
lower energies the dimension should instead run to four, and the dynamics be well- 
described by general relativity. This suggests the intriguing possibility that Nature 
admits a multi-fractal formulation with good ultraviolet (UV) behaviour. 

In the attempt to encode dimensional flow in a physically intuitive framework, in 
[ p3| , [39| , [40H a field theory living in a fractal spacetime was proposed. The key point 



is to replace the standard Lebesgue measure in the action with a Lebesgue-Stieltjes 
measure with anomalous scaling, 

d D x^dg(x), [g] = -Da>-D, (1.1) 

where D is the number of topological dimensions, q is a (possibly very irregular) 
distribution, square brackets denote the engineering dimension in momentum units, 
and < a < 1 is a real parameter interpreted as running with the scale. In the 
ultraviolet, a achieves the critical value 

2 

a ->■ a* := — , (1.2) 

and the measure becomes effectively two-dimensional. Several characteristic features 
were highlighted for the generic measure (|L~T|) [p3| [TO] , but the Lebesgue-Stieltjes for- 



malism is too general to be manipulated. Therefore, the special case of an absolutely 
continuous measure was considered, where dg can be written as 



dg(x) = d xv(x), (1.3) 

for some scalar v. In particular, the Poincare algebra of these models is deformed 
j39fl and at sufficiently small scales an observer should see deviations from Lorentz in- 



variance. At the quantum level, eq. ( |1.2| ) helps field theories with Lebesgue-Stieltjes 
measure to be power-counting renormalizable, gravity included. The reason is that 
the dimension of the couplings changes with the scale and goes to zero in the ultra- 
violet. Dimensional flow also affects cosmology, leading to possible applications to 
inflation and the cosmological constant problem 
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Seen as effective models of quantum gravity, fractal field theories can provide a 
useful tool to describe the physics at scales large enough to smoothen any discrete 
structure into a continuous spacetime, but small enough to retain some properties 
(such as dimensional flow) of the putative full theory. 

Although Hofava-Lifshitz gravity is presently facing some criticism, one of its 
major merits has been to advance the perspective that treating gravity as an ordi- 
nary field theory has nothing wrong in itself, but it requires some rethinking of the 
geometric structure of spacetime. In this spirit, the framework of [^3], |39], [|(J can 
be regarded simply as an alternative to any of the more complicated quantum grav- 
ity scenarios mentioned at the beginning, and the interest shifts from a comparison 
with the continuum limit of discrete models to a detailed study of properties such 
as renormalization and the breaking of ordinary Poincare symmetries. Also, there is 
a conspicuous hiatus between the rigorous mathematical literature on fractals and 
generic ideas of "fractal" structures in field theory and quantum gravity. Letting 
aside dimensional flow, in what sense is spacetime "fractal" ? What is the geometric 
meaning of a fractal in Lorentzian signature PD[? 

It is the purpose of this paper to begin a study of these problems, triggered by the 
qualitative treatment of [39, 40], in a more rigorous framework. Fractional calculus is 
the natural candidate, since it is known to provide a reliable continuum description 
of certain properties of fractals. Fractional measures have the desired characteristics 
of having anomalous scaling and inducing power-counting renormalizability. Their 
study with the techniques of fractal geometry will fulfill the initial expectations and 
unravel a number of new, and perhaps surprising, phenomena. 

The relation between matter and spacetime and the concept of spacetime itself 
are radically different when comparing field theories in Minkowski and general rela- 
tivity. Although realistic models of Nature should include gravity, the introduction 
of the present formalism of fractional dynamics in a flat non-dynamical space will 
help in clarifying how to construct a covariant notion of spacetime with non-integer 
dimension. This will eventually result in a model of multi-fractal geometry with 
non- Euclidean signature (41], j42f . 

We spell out again the physical motivations for studying fractional spaces. We 
can distinguish between direct and indirect applications of these objects to quantum 
gravity approaches. Indirect applications concern the description of certain features 
of quantum gravity such as dimensional flow and the transition from a discrete to 
a continuum spacetime texture. Fractional space is the building block of a more 
sofisticated environment, multi-fractional spacetime [[4~T| [42|]. In turn, the geometry of 
multi-fractional spacetime is an explicit realization of the above-mentioned features. 
The present paper is the starting point of this programme. 

Direct applications of fractional spaces are possible whenever they naturally 
emerge as the limit of some quantum spacetime model. Presently, we are aware only 
of one such example: non-commutative geometry. In P3[, we establish a connection 
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between a certain class of non-commutative geometries and fractional spacetimes. 
The case of K-Minkowski is special, inasmuch as it is obtained as the asymptotic 
regime of multi-fractional geometry at the lowest possible scale. (In multi-fractional 
spacetimes, there exists a hierarchy of scales starting from a fundamental length we 
can identify with the Planck length.) This mapping is realized as an equivalence of 
spacetime measures in both theories. On the non-commutative side, these measures 
preserve the property of cyclicity in the action. Thus, certain non- commutative 
geometries, including the best known one (^-Minkowski), correspond to spacetimes 
with well-defined and precise fractal properties. These spacetimes are nothing but 
fractional spacetimes. Finally, with no reference to other approaches, our programme 
aims at the construction of a UV-finite perturbative quantum field theory. As a field 
theory, one can apply well-known techniques (modulo some modifications) to check 
the desired property of renormalizability. This is illustrated with the example of a 
scalar field theory in [E| fP| . The final goal is to show that perturbative quantum 
gravity is well defined in a multi-fractional setting. 

1.2 Strategy and plan of the paper 

Our programme follows the logic of gradually introducing all the necessary and con- 
cept to formalize the idea of a general multi-fractional spacetime: 

1. First, the fractional analogue of Euclidean space is defined via an embedding 
abstract space (ordinary Euclidean space) and a choice of fractional calculus 
of fixed real order a. In this arena, we can understand elementary notions of 
geometry such as distances, volumes and dimension. 

2. Fractional Euclidean space is a particular model of fractal spacetime. It is 
very important to clarify this concept using tools of fractal geometry which are 
seldom or only partially employed in the literature of quantum gravity. The 
relations between fractional measures and deterministic and random fractals, 
and between different rigorous definitions of dimension, are discussed. The 
Hausdorff (<i H ) and spectral (d$) dimensions of space are non- integer and con- 
stant at all scales. In general, as expected in fractals, 4 ^ 4i although in the 
simplest cases da = ds = Da. 

These two points are the subject of the present paper. There is still little physics 
in this setting for many reasons: the dimension of space is fixed and non-integer 
(in particular, different from 4), the signature is Euclidean, and gravity and matter 
fields are not mentioned. To continue our agenda: 

3. Fractional Euclidean spaces are generalized to fractional Minkowski spacetimes 
with Lorentzian signature. Their symmetry structure is thoroughly analyzed. 
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1. 



A scalar field is introduced as matter, in order to test renormalization properties 
of field theories in fractional spacetimes. 



5. The Hausdorff and spectral dimensions are made scale-dependent via multi- 
fractal geometry techniques. At large scales spacetime is four-dimensional, 
while at small scales it is two-dimensional. Contrary to other approaches 
where phenomenology is the main guiding principle, both these numbers are 
constrained by geometric arguments. 

6. A major advancement in the description of fractal spacetimes is to allow the 
order of fractional operators to be complex. This entails a radical change 
of the physical picture at ultra-microscopic scales, and opens up remarkable 
connections with non-commutative spacetimes and discrete quantum gravity 
models. Discrete symmetries are fundamental, and are progressively averaged 
at larger scales to finally give way to effective continuous symmetries, up to 
the point where ordinary Lorentz invariance is achieved. 

These features were sketched in [f|]J and full details will be presented in the compan- 



ion paper [42||. Finally, as a future work, 



7. One should introduce the notion of fractional manifold in order to realize a 
dynamical fractal spacetime and gravity. Consequences for the cosmology of 
the early universe are expected. 

The strategy, plan, and results of this paper are the following. 

• Section 12. We begin by reviewing the basic tools of fractional calculus in 
one dimension. Our target reader is the physicist unfamiliar with this area of 
mathematics and for this reason we will not attempt a rigorous presentation of 
the latter, preferring intuitive descriptions of its properties. This section is not 
a mere compilation of results taken from the literature. In fact, it also includes 
several proofs of properties of the Caputo derivative, which are rarely found in 
textbooks (section |2.3.2|) , and a novel discussion on mixed operators (section 



2~5l) . Fractional integral measures are a particular case of Lebesgue-Stieltjes 



measure, whose mathematical and physical interpretation is given in section 



2.6. The reader acquainted with fractional calculus can skip this section. 



Section [3|. The first step towards the description of the physical world as 
a field theory on a multi-fractional manifold is to consider a much simpler 
object; namely, a fractional generalization of empty Euclidean space with fixed 
dimension, denoted as £® . This is the ideal playground whereon to understand 
the properties of fractional geometry (section |3.1| ). The notation of fractional 
calculus is not always suitable for this task and it is convenient to introduce 
a new notation (section |3.3|), which makes explicit the geometric properties of 
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E®- After providing hopefully fresh insights on the relations between different 
fractional derivatives (section [3T2D , we define the distance between two points 
(section |3.4j ) . The volume calculations of section [T5] are instructive for several 
reasons. First, they provide examples of multiple fractional integration on non- 
rectangular domains. Second, they highlight how different choices of fractional 
integration affect volume measurements, both in strongly fractional theories 
and in spaces with almost integer dimension. Last but not least, they show 
how the Hausdorff dimension of fractional space is related to the fractional 
order of the measure (section |3.6|) . 



Section We discuss at length in which sense fractional Euclidean space 
is a fractal. On one hand, by considering in detail each of the most charac- 
teristic features of fractal sets, we will show that E® can be regarded as a 
genuine but peculiar fractal, with a continuum structure and anomalous scal- 
ing properties (sections |4.1| - |4.3|) . These are hybrid properties interpolating 



between the trivial case of smooth space M. D (technically a fractal, but with a 
most boring structure) and self-similar sets. By this analysis, we identify the 
symmetries of E® as affine transformations of geometric coordinates, equiv- 
alent non-linear transformations in the embedding coordinates (section |4.2| ). 
On the other hand, fractional calculus can be also seen as an approximation 
of discrete/disconnected fractals under certain assumptions (section |4.4| ). We 
specify the relation between fractional theories and genuine fractal models in 
this section, where the state of the art of this often confusing subject is reviewed 
and clarified. 

• Section [5|. The rigorous definition of the spectral dimension ds on self-similar 
fractals is reviewed. In smooth spaces, it can be found via an operational 
procedure involving a diffusion process. This is the most widely used tool to 
probe "fractal" properties of spacetime in quantum gravity theories, and we 
revisit it at length in section |5.2| . Examining every single ingredient of the 
recipe will allow us to find the spectral dimension of fractional space in section 
f>.3[ As a byproduct, fractional momentum space and the transform thereon, 
replacing the ordinary Fourier transform, are constructed. 

1.3 Original and review material 

Due to the length of the paper, it may be useful to tell apart original from review ma- 
terial. Fractional space E® is a novel object and all related statements are presented 
here for the first time. Many mathematical concepts had been already explored by a 
number of authors, but they had not been collected together into a unified physical 
proposal. We employ the well-established mathematical tools of fractional calculus 
and fractal geometry, and all statements about these two sectors (including defini- 
tions and theorems) are, modulo some exceptions, taken from the literature, which 
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is quoted in the text. Attempting a somewhat artificial discrimination, section |2] is 
a review, except section |2.3.2| (commutation relations are usually reported only for 
the Riemann-Liouville derivative, not for the Caputo derivative) and the text dis- 
cussion in section |2]^. Section [| is original except section |3~T| , eqs. ( p. 1 7|) to ( |3.25|) , 
and the definition ( |3.58| ) of Hausdorff dimension. Section |] is original except general 
definitions and statements about self-similarity (eqs. ([4.i|)-( f4TTl"| ) and ( |4. 18| ) ( P=.2 1| ) ) 
and section O from eq. (4.23) to (4.32). Sections p.l\ and p.2| are reviews, while 



section |5.3| is new, including the discussion of certain conjectures advanced in fractal 
geometry at large. 

Most of previous "fractal" field theoretical proposals assumed spacetime to have 
a non-integer but fixed, non-dynamical dimensionality, thus drawing the attention to 
4 — e dimensions with 0<6<1. On the other hand, here the dimension of space is 
fixed, but in we shall enforce a non-trivial dimensional flow via a superposition 
of fractional measures. Therefore, we shall not be concerned with the almost-integer 
regime except as the infrared limiting case of the multi-fractional theory. A global 
comparison with the early literature on spacetime models in non-integer dimension 
will be done in p2[. 



2. Fractional calculus in one dimension 

We begin by reviewing some basic formulae of fractional calculus. The material in 
this section can be found in the dedicated literature. Among the many excellent 
textbooks on fractional integrals and derivatives, we mention [E5l pj and (471 HH 



(for an historical account, see also ||5]). For reasons which will become clear later, 
here we make extensive use of Caputo derivatives. Multi-dimensional vector calculus 
based on these derivatives has been developed in |50|] (see references therein for works 
adopting other types of integro-differential operators). 

Unfortunately, this branch of mathematics may sometimes produce a feeling 
of estrangement in the reader with field theory/high energy background, thus giv- 
ing her/him the impression of a bizarre topic. This may be due to several rea- 
sons: (i) Contrary to ordinary calculus, there is no unique definition of derivative 
(Riemann-Liouville, Caputo, Erdelyi-Kober, Griinwald-Letnikov, Hadamard, Nishi- 
moto, Riesz, Weyl); (ii) The output of fractional operators is almost always quite 
different from that of ordinary calculus and, in this respect, counterintuitive; (iii) 
Fractional operators are very seldom employed in particle physics, cosmology, and 
physics beyond the Standard Model. 

A deeper inspection, however, shows that: (i) All definitions of fractional deriva- 
tives are related to one another in a precise way (eventually, differences amount both 
to the convergence properties of the functional space on which these operators act 
and to boundary terms in the formulae), so there is no mutual contradiction between 
different operators; (ii) For any choice of fractional operators, the calculus can be 
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defined self-consistently, whatever the look of the formulae; (iii) Fractional operators 
are regularly employed in a number of fields in physics and mathematics, such as 
statistics, diffusing or dissipative processes with residual memory like weather and 
stochastic financial models 0], and system modeling and control in engineering [52 



Moreover, in our context the choice of the Caputo operators over the others will be 
motivated by independent arguments, thus further neutralizing issue (i). 

Although by now some particular symbols are widely employed, for the sake of 
simplicity, and in order not to overburden the unfamiliar reader, we will reserve a 
simplified notation for the operators we use the most. 

2.1 Integrals and derivatives 

Let 

x G [x ,xi\ (2.1) 

be a real coordinate variable defined on an interval with constant extrema xq and 
Xi, which may be taken to infinity if desired. We define a space of functions f(x) 
on this interval, such that all the following integro-differential operators will be well 
defined. A space of particular interest is AC u [xq, x\], that of functions which are 
absolutely continuous on [x , Xi] up to their n — 1 derivative. This is equivalent to the 
space of Lebesgue-summable functions with summable derivatives &*, j — 0, . . . , n, 
almost everywhere in the interval. Another space we shall use more or less implicitly 
is L p (xq,Xi), that of Lebesgue-measurable functions / on [xo,xi] with finite p-norm 
From now on, we assume that for every mathematical statement the functional 



space is suitably chosen. One can find precise theorem declarations in [[|5|, |4q[ . 
Let / G L p (xq,Xi) and let 6 be the Heaviside distribution: 

6(x) = { 1, X> ° . (2.2) 



,0, a;<0 

We introduce the left fractional integral of order a as 



(2.3) 



Here a G C is a complex constant parameter, which we shall restrict to be real for 
our purposes. This formula is naturally suggested as a generalization to non-integer 
n of the Cauchy formula for the n-times repeated integration: 

rx ryi rVn-i 

(I n f)(x) = dy x dy 2 -- dy n f(y n ) 

J XQ J XQ J XQ 



1 



n-11! /„ 



dx' (x - Z 7 ) /OO • 

XQ 
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One can also define the right fractional integral of order a 

dx' 



(2.4) 



where integration now is from x to the end of the interval. 1 Because of the x depen- 
dence in the step function, the output of fractional integrals is a function of x. These 
operators are bounded if / G L p (xa,xi), for every 1 < p < +00 [46], eq. (2.1.23)]. 
Sometimes we will need to specify the integration domain explicitly and we shall 
denote the left and right integrals with I a = I x x and I a = I XX1 - 

In parallel, for any / G .AC n [xo, #i] the left and right Caputo derivatives of order 
a ]53], exist almost everywhere in [a?o,Xi] 0, theorem 2.1]: 

(d a f){x) := {I n - a d n f){x), n-l<a<n, (2.5) 
1 f Xl dx' 

r(n - a) J xo (x - x') a+1 -« V ' x K ' K J 

(d a f)(x) := (I n - a d n f)(x) , n-l<a<n, (2.7) 

1 ■ -ev-x)d$fw, (2.8) 



T(n - a) J XQ (x' - x) a+1 - n 

where d is the ordinary first-order partial derivative and n > 1 is a natural number. 2 
We shall be interested in the particular case 

< a< 1 (n = 1) , (2.9) 

which simplifies the above expression for the left derivative as 
1 f Xl dr' 

a similar statement holds for the right derivative. An elementary but typically ig- 
nored fact is that, under the transformation 

X — ¥ Xq + X\ — x , (2-11) 

left operators are mapped into right operators: 

(I a f)(x) = (I a F)(x + Xl -x), (B a f)(x) = (d a F)(x + Xl -x), (2.12) 



1 The left fractional integral is variedly indicated in the literature with the symbols Ko I", Xo D x a , 
x d~ a , The right integral, with x /°, X D~^ x d~?, 

2 The left fractional derivative is variedly indicated in the literature with the symbols Xo D", Xo dxi 
D XQ _. The right derivative, with X D" ± , x d Xl , D% 1+ . All these symbols are further decorated with 
tilde's, in bold font, or with the superscript C when they denote the Caputo derivative. Fractional 
operators are often called differintegrals, since one can analytically continue derivative expressions 
to a < and vice versa. 
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where F(x) := f(xo + X\ — x). Therefore, it is sufficient to study the properties of, 
say, left operators, and infer their right counterparts by using eq. Q2.12 ). For this 
reason, in the following we mainly concentrate on left operators. 

An alternative, inequivalent definition of fractional derivation is obtained by 
exchanging the order of integration and derivation. This corresponds to the left and 
right Riemann-Liouville derivatives (n — 1 < a < n) 

Ud a f)(x) := (d n I n - a f)(x) 

1 C X1 dr' 

d n x t ^——e(x-x>)f(x>), (2.13) 



T(n -a) x J XQ (x - x') a+1 - n 



L B a f)(x) := (d n I n - a f)(x) 



' d n x f * - d *'_J (x' - x)f(x') . (2.14) 



r(n - a) x J XQ (x' - x) a+1 - n 

There is a precise relation between Caputo and Riemann-Liouville derivatives. If 
d J x f is continuous on [x ,x] for j = 1, . . . , n, one has 

n-l 

(d a f)(x) = Ud a f)(x) - T p"° r (&f)(x ) , n - 1 < a < n. (2.15) 



(x — Xo) J a 

In particular, for n 



(d a f)(x) = Ud a f)(x) - {X Xo) ~ /(s„) , < a < 1 , (2.16) 

1 (1 — a) 

and the two derivatives are the same if f(xo) = 0. In general, the two types of 
derivatives differ in the boundary conditions. While {d a f)(xo) = by definition, 
because of eq. ( |2.15| ) 

Ud a f)(x ) = & (d>f)(x ) = 0, j = 0,...,n-l. (2.17) 

So, the Riemann-Liouville derivative of a constant is not zero. Right derivatives 
obey similar relations. 

2.2 Examples of fractional derivatives 

Let us see some examples of calculation with Caputo left fractional derivatives with 
n — 1 < a < n. A simple function to consider is a power law, 

f(x) = (x - x*f , x*e[x ,xi], (3eR. (2.18) 

If(3 = m = 0,l,...,n — 1, then both the left and right Caputo derivatives vanish for 
every x*, 

d a (x -x*) m = = B a (x -x,) m , m = 0,1,..., n-l. (2.19) 
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Barring these cases, and assuming xo ^ — oo, one has 



r(/3 + i) f J( (x' - x*)^-« 

' ax 



T(n - a)Y{fi - n + 1) J Xo (x - x') a+1 ~ n 



r(n - ot)r(/3 - n + 1) 7o [0 - s ) - y]^ 1 -" 



r(/? + i; 



r(n - a + l)r(/3 - n 



(x - X*) 13 n (x - Xq) 



x 2 F 1 (l,n-/3;n+l-a,- — —J , (2.20) 

\ X* — Xq J 

where we used pH eq. 3.196.1] and 2 -Pi is the hypergeometric function. In the special 
case x* = xo, we get eq. Q2.19| ) and [^5], eq. 3.191.1] 

d a (x - xo)' 3 = -I^±j- (x - x )^ , /? ^ 0, 1, . . . , n - 1 . (2.21) 
1 (p — a + 1) 

When the lower extremum is Xo = —00, one can employ eq. 3.196.2 of 



™Aa+l— n 



r(n — a)T(/3 — n + 1) y_ 00 (x — x' 

;-i)^ w r(/3 + i) /-+ 00 (y + x,)' 3 -" 

ay 



r(n — a)T(/3 — n + 1) (y + x 

(-ipr^ + i)r( a -/3) 



r(l - n + 0)T(n - p) 

F(p + 1 — a) sm[7r(p — a)\ 

where we used the property T(z)T(l — z) = tt/ sin(7T2;). This expression is ill defined 
for (3 = a. Otherwise, it is real under certain conditions on the values of a and 
and the sign of x — x*; consistently, it vanishes for (3 = m. The Riemann-Liouville, 
Caputo and Griinwald-Letnikov derivatives all collapse to the same operator when 
Xo = —00, the Liouville fractional derivative ^d" (in particular, oo<9 a l = 0). When 
regarded as an approximation in the limit t ^> t , this operator is employed in 
mechanics to describe "steady state" systems, that is, systems which evolved well 
after the initial transient phase at time to- The right derivative (or integral) with 
Xi = +00 is called Weyl derivative (or integral). To get the Weyl differintegral from 
the Liouville differintegral, it is sufficient to set xo = —x\ in ( |2.12| ) and then take the 
limit xo — > —00: 

U-PiOO) = Ul a F)(-x) , Ud a f)(x) = (^FX-x) , (2.23) 
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where F(x) := f(—x). 

The a-th order integral of (x — x*Y is given by the analytic continuation a — > —a 
of the above formulae for any /3. In particular, to prove some commutation theorems 
we will need the integral eq. (2.1.16)] 



I a (x — Xo) 



r(/3 + i 



T(f3 + a + l) 
The extension of eq. ( p.22 ) to negative a is 



(2.24) 



r a. 



dx 



, (x' - x^Y 



r.l\l-a 



y=-x 



-1) 



oo 

13 r+oo 



r a 



dy 



x — x \ 
(y + a;) 1 "" 



r(-/3) 

4 )-r(-a-/3) 



(x - x*f +a 



Another example is the derivative of the exponential function, 

f( x ) = e Xx J O^AGl. 

For finite Xq, we have 



(2.25) 



(2.26) 



a „\x 



d a e 



A' 



Tin — a) 



y=x -xo 



^n^XxQ rx—xo 



e \x 

, ^ ix - x') a+1 - n 



Yin — a) 
Tin — a) 



dy 



[(x - x ) - y} a + 



1-n 



e X/ y[n — a, X(x — x )] 



(2.27) 



where we used [55], eq. 3.382.1] and 7 is the incomplete gamma function. Comparing 
its series representation 



7(6-1,2) = r(6- l)z b - l e- z Y^ 



fc=0 



T(k + b) 



with that of the two-parameter Mittag-Leffler function [q6 



E a , b (z) := J2 



k=0 



T(ak + b) 



a > 



one can also express eq. Q2.27| ) as 

d a e \x = X "e Xxo (x - x ) n - a E ltn+1 _ a [X{x - x )} 



(2.28) 



(2.29) 



13 



Since j(z, +00) = T(z), the exponential is an eigenfunction of the Liouville 
derivative (xq = —00), with eigenvalues X a : 



This can be obtained also from [55, eq. 3.382.2]. For finite xq, the eigenfunction of 



the Caputo derivative is the one-parameter Mittag-Leffier function E a (z) := E a ^(z)\ 

d a E a [X(x - x ) a ] = XE a [X(x - x ) a ] , (2.31) 

which stems from differentiating eq. ( j2.28| ) term by term via eq. (|2.21[) (and remem- 
bering that the constant term k = gives zero). Thus, the Mittag-Leffier function 
E a [X(x — Xq) 01 } can be considered as the fractional generalization of the exponential 
e z = Ei(z). With the same procedure, one can get other expressions such as 

d a {(x - x fE a ^ +1 [X(x - x ) a ]} = X(x - x f- a E a ^ a+1 [X(x - x ) a ] . (2.32) 

Notice, however, that one should exercise care in inferring some results by analytic 
continuation of others. For instance, one does not recover eq. (|2.31|) from ( |2.32| ) in 
the limit a — > a, /3 — > 0, because E a>1 _ a (z) = [T(l — a)]~ l + zE a (z). The first term 
would give an extra contribution (x — £o)~ a /T(l — a) to eq. ( |2.31| ). The reason is 
that, contrary to the case of the Riemann-Liouville derivative, the limit — > of 
eq. ( |2.21| ) does not give the correct result ( |2.19| ). In fact, when differentiating the 



left-hand side of (|2.32|) term by term one picks up also the k = contribution, which 
is zero if /3 is set to zero from the beginning. 

Analogous formulae can be obtained for the right derivative by making use of 
eqs. (|2.12| ) and ( |2.23|) . For example, the right version of eqs. Q2.21|) , (|2.22|) and ( P.3Q ) 



are simply 

d a { Xl -xf = r ff + 1 } -A Xl -xf-\ /3^0,l,...,n-l, (2.33) 
1 (p — a + I) 

T(p + 1 — a) sm[7r(p — a)\ 
2.3 Properties of fractional operators 

We list some properties of fractional operators, focussing on left Caputo differ- 
integrals but quoting some results also for Riemann-Liouville operators. Here, 
n — 1 < a < n and m — 1 < < m. 

We detail the derivation of most of the results for two reasons. One is that some 
of the theorems below on Caputo derivatives cannot be found in [ 45 , [16|] . Another is 



to give the unacquainted reader a few examples of the subtleties of fractional calculus. 
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In order to prove some of the statements, we will invoke the first and second theorems 
of fundamental calculus on the interval [xo,xi]: 

(d n I n f){x) = f(x), neN. (2.36) 

n— 1 , 

(I n d n f)(x) = f(x) - ~ (x - xoYi&flixo) . (2.37) 

2.3.1 Limit to ordinary calculus and linearity 

When a = n G N, one recovers ordinary calculus of integer order n |55|, eq. (2.4.14)]: 

lim d a = d n , lim B a = (-l) n d n . (2.38) 

a—*n a—tn 

In particular, when a = derivatives and integrals collapse to the identity operator. 
Fractional operators are linear: 

O a [c 1 f{x) + c 2 g{x)} = c 1 {O a f)(x) + c 2 {O a g){x), O a = d a ,B a ,I a ,I a . (2.39) 

2.3.2 Commutation relations 

The commutation relations 

QolqP = QPQ* ^ O a>0 = Qa ^ Qa <0 = ja ? ^.40) 

are valid under one of the following conditions. 



jaj/3 _L jPja 

If a, (3 < 0, at almost every point x G [x ,Xi], fractional integrals obey the 
semi-group property (J45], eq. (2.100)]; []46], eq. (2.1.30)]) 



jajfi = jfija = ja+f)^ V tt, /3 > . (2.41) 

This property holds everywhere in [xo,Xi] if a + /3 > 1. 
U/3>a, 

= I?- a , (2.42) 
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while if (3 < a we have, letting k — 1 < a — (3 < k (so that k = n — m < n), 



" la; Xoj (a*/)(x ) 

— i I .1 ■■/■■-.> i 

J = 



1 (**-'/)(*) + E r(i + TlT-a) {djf){xo) • (2 ' 43) 



J 

On the other hand, when (3 > a, 



(I?d a f)(x) H (I r- a d n f)(x) 

= i^- a rd n f{x) 



m - 1 ra;r+7-l) ( ^ /)(Xo) ■ (2 ' 44) 



while for /3 < a 

(I^d a f)(x) = {I (3 r- a d n f){x) 

(d a -^I n d n f){x) 

n-l 



J* (d^f)(x) - J2 J* *°r«~ a S (&f)(xo) . (2.45) 

~ r(j + i + p - a) 

Therefore, eq. ( |2.45| ) is valid for all a, (3 > and 

d a I p = I^d a , a,(3>0 (d j f)(x o ) = 0, j = 0, . . . ,n - 1 . 

(2.46) 

In comparison, for a, /3 > the Riemann-Liouville relations are ^d ! 13 = d a ~^ 
(for any a and /3, not just a > (3) and 

(I^f)(x) = ( RL d^f)(x) - £^y^yM a -V)(x ) , 

and, via eq. ( PT7|) , one has RL d a I^ = l\ L d Q if (d j f)(x ) = 0, j = 0, 1, . . . , n-l. 
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Qa Qm _i_ Qm Qa 

If = m E N+ holds if 

(^/)M = o, 

In fact, on one hand we have 



j — n, n + 1, . . . ,n + m — 1 



(2.47) 



Q a Q m — jn—ocQnQm 

j(n+m) — (a+m) Qn+m 
_ Qa+m 



(2.48) 



On the other hand, by using eq. ( |2.15| ) twice and d m RL d a = RL d a+m , one gets 

n+m— 1 



(d m d a f)(x) = (d a+m f)(x) + J2 YJi 



X — Xq 



tj—m—a 



]=n 



+ j — m — a 



i&f){x ), (2.49) 



hence the result. For the Riemann-Liouville counterpart of eq. Q2.40D , plug 
(pl5|) into ( gggD : then, 



m— 1 



( RL ^«9"7)(x) = ( RL <9«+"7)(x) - E fTl 



jj— m— a 



3=0 



(1 + j — m — a) 



(37) (*o) 



and RL d a d m = d m RL d a when (d^ f)(x ) for j = 0, . . . , m — 1. 

Another case of interest is the commutation relation between two Caputo 
derivatives. When a, (3 > 0, one has 



(d a d^f)(x) (I n - a d n d?f)(x) 



n+m— 1 



jn—at T „ y— n- 

(/ n- Q ^ + n /)(x)+ £ _ ^ *0j _ (ffl - /)(xo) 



2.24) 



j=m 
n+m— 1 



Xl+J-n-/3) 



x — x 



(SVXso) 



x - xo)^^ 



i=o 



(l+j-a-/3) 



(^7)(x ). (2.50) 



(Incidentally, notice that one cannot analytically continue eq. ( p. 50 ) to the 
cases P — m and a = n. For = m, there is an obstruction in the step from 
the first to the second line, while for a = n there is an obstruction from the 
third to the fourth line. The correct expressions ( |2.48| ) and ( [2.49 ) are obtained 
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from the first and third line, respectively.) Switching (at,n) and (/3,m) and 
comparing the expressions, one finds that 

gagP = gPg* 

if either 

n = m, (2.52) 

or 

a + p = n + m-l, (2.53) 

or 

(x ) = , j = r , r + 1, . . . , f - 1 , (2.54) 

where r = min(n, m) and f = max(n,m). Interestingly, fractional derivatives 
do commute if a and (3 have same integer part. The reason is that eq. ( |2.50|) 
is symmetric in a and 0. Notice also that d a d /3 ^ d a+ ? unless (&>f)(x ) = 



for j = 0, . . . , r — 1 or eq. (|2.53 ) holds. In the latter case there fall the values 
a — j3 = 1/2, and one has d^d^ = d. 

The commutation relation for the Riemann-Liouville derivatives is ( [[45| eq. 
(2.126)]; H, eq. (2.1.42)]) 



(x — X ) j a 



jr[ r(i -3 -a) 

Switching (a,n) and (/?, m), the relation 

f)a gP _ a/3 aa 
Rh u «L U — RL U RL U 

holds if, simultaneously, (rlCT^/X^o) = for j — 1, . . . , m and (rlC^-' /Xxo) 
for j = 1, . . . , n . Using eq. ( |2.17| ), the combined condition is 



(07) (*o) = 


0, j = 0,l,...,f-l 


a, /3 < 1 and m 


G N+, 






Qaj/3 = jfSga 


& f(x ) = 0, 


d a d m = d m d a 


(df)(x )=0, 


gagf) = gPg* ; 


V a, (3 



(2.55a) 
(2.55b) 
(2.55c) 
(2.55d) 

for the Caputo left derivative, while 

^d a lP = l\ L d a & /(x ) = 0, (2.56a) 

gagra = gm ^ga ^ = q ? (2.56b) 

,d a RL dP = RL d? RL d a & f(x ) = (2.56c) 



RL ( - 



R L L 
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for the Riemann-Liouville left derivative. Finally, the Liouville differintegral oo<9 a 
always commutes, for sufficiently good functions / (i.e., continuous with continuous 
derivatives and which fall to zero with their derivatives sufficiently fast for xq — > — oo) 

[0: 



r 3 



d a I 13 



t/3 ja 
-J 1 oo- 1 

TP f)<* 
oo- 1 oo u 



I J- 



d 



,a-/3 



OO u OO*-' OO u OQ*-S OO ^ 



V a, /3 > 0, 

V a, (3 > 0, 

V a, (3 > 0. 



(2.57a) 
(2.57b) 
(2.57c) 



2.3.3 Fundamental theorems of calculus 

The fractional derivative is the left inverse of the integral. Setting a = (3 in eq. ( 2.42 ), 

(d a I a f)(x) = f(x), a>0. (2.58) 

This equation holds also for the Riemann-Liouville derivatives (|45|, eq. (2.106)]; [f46| , 
eq. (2.1.31)]). 

The extra term in eq. (|2.15|) is responsible for the following, important difference 
between Riemann-Liouville and Caputo derivatives. One of the points where frac- 
tional calculus may show its worst trickiness is upon generalization of the theorems 
of calculus, such as the Newton-Leibniz formula 



dx'(df)(x') = f(x)-f(x ) 



•ro 



The same formula is not valid for the Riemann-Liouville derivative. In fact, one can 



show that (P, eq. (2.108)]; g§, eq. (2.1.39)]) 



a-j 



id n ^r- a f)(x ) 



n — 1 < a < n , 



so that for n = 1 

(I a RL d a f)(x) = f(x) 



T(a) 



< a < 1 . 



On the other hand, the Caputo derivative is the only fractional derivative obeying 
a simple Newton-Leibniz formula without imposing particular boundary conditions 
on the functional space. Thanks to eq. ( 2.41 ), one has 

(I a d a f)(x) = (I a I n - a d n f)(x) = (I n d n f)(x) , 

so that we get 



n— 1 



(i a d a f)(x) = f( x ) -J2^(x- x y(&f)(x ) , 



n — 1 < a < n . 



(2.59) 
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which could have been obtained from eq, ( |2.45| ) with a = (3. For n — 1, 

(I a d a f)(x) = f(x) - f(x ) , < a < 1 . (2.60) 

Therefore, both fundamental theorems of calculus are satisfied by the Caputo deriva- 
tive. If the terminal points are not equal in I a and d a , composition laws become 
more complicated; we do not consider this case, since these operators are always 
thought of as defined on the same domain. 

2.3.4 Leibniz rule, non-locality, composite functions 

Unfortunately, the Leibniz rule of derivation for a product of functions / and g is 
complicated whatever the choice of derivative, since it contains an infinite number 
of terms g5|, eq. (2.202)]. If f,g e C°° in [x Q ,x], then 

jw-tQvf*^,). c a Ka-y + a + D - ™ 

where <9 Q_J = P~ a are actually integrations for j > 1. If / and g are, respectively, 
analytic and continuous in [xq, x], eq. ( |2.61| ) is valid also for a < (i.e., for fractional 
derivatives replaced by fractional integrals) and for Liouville/Weyl operators 



Equation ( p.61|) shows the non-local nature of fractional operators: fractional 
integration by parts or derivation gives rise to an infinite number of terms. Setting 
g = 1 or g = 9(x — Xq), one obtains an expression of the Caputo derivative d a as an 
infinite series of ordinary derivatives. For n — 1, 

1 f{x) - f(x ) 



(d a f)(x) 



r(l-a) (x-x ) a 



where we used eqs. ( ^.16[ ) and ( |2.24| ). 



Equation ( ^.62| ) is useful for writing down the fractional derivative (or integral) 
of a composite function f[g(x)]. In fact, the simple rule 

W(9)] = j- 9 *9 (2 ' 63) 

no longer holds. The j-th derivative d ] x [f(g)] can be further expanded as a series, 
using the Arbogast-Faa di Bruno formula [||5|, section 2.7.3]. 

2.3.5 Integration by parts 

For an ordinary integral where the integrand contains fractional derivatives, one has 



H eq. (2.1.50)] 



d x fd a g= / dxgd a f, (2.64) 

XQ J XQ 
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and the same formula, under suitable conditions on the functions, holds for d a and 
d a replaced by I a and I a , respectively. One can also consider the opposite situation, 
i.e., a fractional integral where the integrand has only ordinary derivatives: 



1 f Xl dx 



X'l 



dx Gdf G 



r(a) (xi -xy- a 



= - dx fdG 


J XO 






rxx 


r(a) J 


xo ( X l 


= -r \ 

x ,x 1 1 


fdg + 



dx ( 1 - a 
;f[dg + g 

X\ — X 



\X—a • 



fg . (2.65) 



The integration by parts of fractional integrals of fractional integrands can be inferred 
by combining these two cases. As an example, we take a fractional integral over the 
positive semi-axis. Denoting with a left subscript the lower (upper) extremum of 
integration in left (right) derivatives, 



I r+oo 

Io,oo {gd P f} = ^r-J dxx a ~ l g^f 
r(a) Jo 



63) 1 



+oo 



dxf(x) 00 B^[x a - 1 g(x)}, (2.66) 



o 



I» 

assuming that the functions /, g are good enough for all steps to be well defined. 
2.4 Exterior derivative 

Fractional differentials were early proposed in [57|-[S0| for the Nishimoto derivative, 



in 



61 ] [63] for the Riemann-Liouville derivative, and in [^, 64, 65] with Caputo 
derivatives; early applications of fractional differential forms to mechanical systems 



can be found in |65|-||69||. We shall mainly follow the results of [5(| |61|, with adap- 
tations. 

The fractional exterior derivative is defined via the left derivative as 

d a ■= (dx) a d a , (2.67) 

and a right definition also exists: d a := (dx) a B a . The reader may wonder, on one 
hand, whether this definition is natural and, on the other hand, about the meaning 
of the writing (dx) a , the "a-th power of dx." These two questions are actually 
interrelated. We shall postpone the answer to the first in section |3.3| , where we will 
see how the second fundamental theorem of fractional calculus ( |2.60| ) immediately 
suggests eq. (j2.67|) as the obvious candidate for the fractional differential of a function. 
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As for the second question, the object (dx) a is a compact rewriting of the fractional 
differential of a certain function of x, which will be later recognized as the natural 
coordinate in fractional space. By eq. ( |2.21| ), one sees that d a (x — x ) a = T(l + a) 
and 

d a (x -xo) a = T(l + a){dx) a 1 (2.68) 
and (dx) a is shown to be, up to a constant, the fractional differential of 



Q 



r(i + «) 



-a . 



(2.69) 



Therefore, we can recast eq. ( |2.67| ) as 

d a = d a qd a . 



(2.70) 



For integer a, the fractional differential behaves as the ordinary one. Taking again 
eq. ( |2.21|) , one has d°(x— x^Y — { x ~ i d 1 {x—xo) 13 = (3(x—xo) ,3 ~ 1 dx, d 2 (x—xo)^ = 
(3((3 — l)(x — x ) l3 ~ 2 (dx) 2 , and so on. 

Notice that the left-hand side of eq. ( |2.68| ) seems to be ill defined for x = Xq or 
Xq = —oo at any x, but the right-hand side just shows that these are artifacts of the 
presentation in fractional coordinates. The internal structure of d a conspires with 
that of the fractional coordinate to give a finite, well-defined result. 

An exact fractional 1-form is the differential of a scalar function, d a f. A generic 
fractional 1-form is u> = (dx) a f(x), for some function /. The exterior derivative 
of uj yields a 2-form, which requires an extension of the coordinate space to many 
dimensions. In that context we shall describe a geometric interpretation of fractional 



differentials [^9 



2.5 Mixed operators 

After choosing to work with the left or right sector, one can consistently define all 
the elements of differential calculus within the same sector. However, in a fractional 
calculus of variations we expect to have a mixing of the sectors because of eq. (|2.64|) . 
This may be a first reason to also consider versions of fractional calculus with mixed 
operators. 

As far as derivatives are concerned, Cresson defined the complex linear combi- 
nation JT(J 3 

iA + 1 



-d a + 



iX 



-d a 



such that integration by parts becomes 



1 dxfV$g = - I 1 dxgV a _J. 

XQ Jxo 



(2.71) 



(2.72) 



3 This operator was actually denned for Riemann-Liouville derivatives of mixed fractional order. 
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When A = — i and A = +i, one recovers d a and — d a , respectively. For A = 0, the 
same derivative operator would appear on both sides of the equation. However, the 
operators d a and B a have complementary domains and it is not possible to define 
a generalized fractional coordinate associated with the operator V". In the interval 
[xo,x), the natural fractional coordinate is eq. ( |2.69| ), while in the interval (x,xi\, it 
is 

(xi — x) a 

« '■= n I \ ■ 2 ' 73 
1 (1 + a) 

One of the extrema in differintegral operators varies and, very roughly speaking, one 
cannot envisage a functional which is constant in the first and second interval for 
any x (compare eqs. ( 2.21 ) and ( p. 33 ) with = a). There seems to be no natural 



geometric interpretation of a theory defined with T>\, for any A. 

Another motivation to construct mixed operators is in the way fractional inte- 
gration is carried out. In fact, there are different prescriptions for generalizing the 
definite Lebesgue integral 

dx'f(x') 



J x 



to a fractional integral. The left operator J Q , eq. ( p.3|) with x = x±, carries a measure 
weight (xi — x')"" 1 , such that the major contribution to the integral comes from 
the area under the curve with x ~ x\. On the other hand, the right operator I a , 
eq. (|2.4 ) with x = xq, carries a measure weight (x' — xq) 01 ^ 1 dominating near the 



lower extremum, x ~ Xo: the output value of the two integrals is not the same. This 
will not result in different physical qualitative properties, but will eventually lead to 
different physical measurements, as we shall see in the next section. Therefore, it 
may be interesting to explore other possibilities. A mixed-type generalization of the 
definite integral draws inspiration from the splitting 

x\ rx* rxi 

dx'f(x')= / dx'f(x') + / dx'f(x') , V x* G [xo, xi] , 

' XQ J XQ J X* 

so that one can define 

Itf:=[(I a + I a )f}(^), (2-74) 

or, more explicitly, 

^XQ,x t ,xif ■ ^XQ,x t f Ix*,xif • (2-75) 

Actually, and contrary to the Lebesgue case, due to the fractional weights the split- 
ting is no longer arbitrary and J™ is a class of inequivalent integrals parametrized 
by x* e [x ,xij. 

Notably, I a can be written in terms of just one sector [p5|, [7T], After 
obvious coordinate transformations, one has 

I rx*-xo a t i i rxi-x* » 
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Upon changing integration variable in the first term as x' = [(x* — xq)/(xi — 
we obtain 



T a f = T a f 



where 



Thus, 



I a f 



f+(x) ■= f(x* + x) 



ja p 

0,xi—x*J - 



T(a) 



f-( x ) :=/(&*- X * - X ° x 



(2.76) 
(2.77) 

(2.78) 



where 



U( x ) + 



(2.79) 



(2.80) 



Therefore, the mixed integral I a is equivalent to a right integral on the positive semi- 
axis (the weight dominating near the origin), acting on a modified function space. 
Expression ( |2.79| ) simplifies under some conditions. First, if the definite integral is 
symmetric (xq = —x\ = —R), we have 



I-R,x t ,Ftf 



r(a) 
1 



x 1 



f(x* + x) + 



R + x* 



,-. * R ~\~ 



(2.81a) 
(2.81b) 



In particular, the symmetric choice x* = yields 

2 



* f 

-R,0,RJ 



1 



x 1 



im + /ho] 



(2.82a) 
(2.82b) 



When / is even, f(x) = f(x). We will use this property later to calculate multiple 
volume integrals. Also, sending R — > +oo in eq. Q2.81 ), the mixed fractional integral 
can be presented as 



Im.,x t f 



T(a) 



f( x ) = 2 \-f( x * + x ) + f( x * - x ) 



(2.83a) 
(2.83b) 
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which can be further specialized to x* = 0, if desired. Notice that this is a gener- 
alization of the ordinary integral over the whole real axis, despite the fact that the 
presentation ( |2.83| ) is on the positive semi-axis. Taking the absolute value \x\ in the 
measure, one can formally extend this presentation to the whole axis, but we prefer 
to keep ( |2.83| ) because it makes explicit the existence of a boundary (the special point 
x = 0). This will be the object of much discussion when defining the symmetries of 
fractional spacetime W2 . 



Unfortunately, also the mixed integral has unattractive properties. Integration 



by parts follows from the results of section [2.3.5| . When the integrand has only 
ordinary derivatives, one has 



II {gdf} © -r^\fdg + ^—^fg 



(-1) - 1 , f(x)g(x) 

+ v ' . — lim / v m / . 2.84 
r(a) (x* - x) l - a v ' 

The last term must vanish, thus imposing a function space such that for any / and 
g, (fg)(x) ~ (x* — x) 1_cH_e near x ~ x*, where e > 0. The physical meaning of this 
constraint is not clear a priori. Finally, the fundamental theorems of calculus do not 
hold, since there is no left inverse of the mixed integral (a linear combination of left 
and right derivatives would produce cross-terms which, in general, do not cancel). 

2.6 Lebesgue— Stieltjes measure and interpretation of fractional integrals 

Equation ( |2.69|) can be also regarded as the Lebesgue-Stieltjes measure associ- 
ated with left fractional integrals. Indeed, the distribution q(x) defines a mea- 
sure Q a over the interval [xo,x]. This measure is the Caratheodory extension of 
m e ((x*,x]) := g(x) — g(x*), m g ({xo}) : = 0, for any x* 6 [xo,x]. In fact, g a is mono- 
tonic, non-decreasing and right-continuous. Furthermore, the properties of measures 
are satisfied. We can see this intuitively by considering m g : (i) m^(0) = 0, (ii) 
^Iq(A) < m g (B) if A G B ^ [xo,x], (iii) given a countable or finite union of sets, 
m Q {\J i A, l ) < J2i m e(Ai), where the equality holds if the Ai are disjoint Borel sets; in 
particular, m Q (A\B) = m Q (A) —m g (B). Therefore, the Riemann-Liouville integral 
( p.3|) can be regarded as a Lebesgue-Stieltjes (or Radon) integral [[73]— [[75]: 



r*f = P dg a (x') /(*') , g a (x') = ~ %r^ , (2.85) 

where we made a slight abuse of notation and identified Q a (x) with m e ((x , x\) = q(x). 
An important feature of the measure is the scaling property, inherited from q, 

Qa (Xx) = \ a g a {x) , A>0; (2.86) 

namely, the measure of the set obtained by a rescaling x — > Xx is the same as the 
original set, multiplied by a factor A". 
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Fractional integrals admit neat geometrical [76, 77] and physical [[7g, 79 inter- 
pretations. Consider a function f(t) and the time integral 



(I a f)(ti 



to 



h dt^-^f(t) 

Via) 



to 



dg a (t)f(t). (2.87) 



to 



The geometric meaning of the left fractional integral (|2.87|) with a ^ 1 fixed is 
shown in figure [I]. The continuous curve in the box is given parametrically by the 
set of points C = {(£, g a {t), f{t))}, where / is some smooth function. Projection 
of C onto the t-f plane (g a = const.) gives fit), while projection onto the t-g a 
plane (/ = const.) yields g a (t). Now, build a vertical "fence" under the curve C, 
and project it onto both planes. On the t-f plane, the shadow of the fence is the 
ordinary integral, 



dtf(t) 



(2- 



to 



On the g a -f plane, the shadow corresponds to the fractional integral Q2.87|) , the area 
under the projection of C on such plane. 



/if) 




Figure 1: Geometric interpretation of Lebesgue-Stieltjes integrals as "shadows" of a 
"fence." fit) is a generic smooth function and the measure g a (t) is given in eq. ( 2.85| ) 
(specifically, in the figure a = 1/2). 

The behaviour of the measure weight v a (t\ — t) leads to different physical sce- 
narios at the extreme values of the interval < a < 1. We can regard this weight 
as a memory function and the fractional integral as a memory flux. If a — 1, the 
fence lies in the t-f plane and the memory function v\ = 1 equally weighs all the 
points from the initial time to to the final time t, eq. (|2.88|). Processes described 
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by integer integrals retain all the memory of the past history. The integral has the 
usual meaning of "area under the curve fit) in the interval [to,ii]." Also, if fit) is 
the speed of a point particle, I 1 f is the operational definition of the distance covered 
in the time interval At = t\ — t - 

When a = 0, the limit of the memory function in the sense of distributions is a 
delta, v a ^o(ti —t) — 5(ti — t), the fractional integral becomes the identity operator, 
and the integral of a function from t to t\ equals the function itself, evaluated at 
the final point t\\ 

Yimil a f){t l ) = f{t l ). (2.89) 

The past history is completely forgotten. Systems with no memory retention are 
called Markovian, and are well described by fractional calculus in the limit a — > 0. 
Thinking of information as carried by "states," a roughly corresponds to the fraction 
of states preserved at a given time t. In turn, loss of information corresponds to a loss 
of energy, at a rate 1 — a. Hence, fractional systems are dissipative jBlfl . Examples 
are percolation clusters, porous media, collision systems, and Brownian motion. The 
one-sidedness of fractional operators, in fact, is responsible for the irreversibility of 
time fS|. 



3. Fractional Euclidean space 



The extension to D topological dimensions (i.e., to D coordinates, where D G N + ) 
is straightforward. Each direction is associated with a fractional "charge" a M . The 
corresponding Lebesgue-Stieltjes measure is 



D 

Qa{x) = (^) Qa^) , 



(3.1) 



with rescaling 



o a (Xx) = X^^Qaix) , A>0. 



(3.2) 



In general, the D parameters a M can be different from one another. To make the 
presentation as simple as possible, we shall make an "isotropy" assumption: namely, 



a. 



a . 



V/i = l, 



D 



(3.3) 



This restricts the analysis to fractional manifolds (section |3.3|) where the fractional 
charge is equally distributed among the directions, and all of them are treated on 
an equal footing. Anisotropic configurations are possible and different choices of 
the set {a^} (modulo permutations) correspond to inequivalent geometries. Some 
anisotropic measures were given in [j39|. 
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Before discussing fractional Euclidean space, we need the extension of fractional 
differentials to many dimensions. The partial fractional derivative along the /i direc- 
tion is 

3?:=^, [5?] = "- (3-4) 
A simple summation over coordinates yields 

d a ■= (dx^d* , [d a ] = . (3.5) 

The Einstein convention of summing over repeated upper-lower indices is employed. 
Arbitrary fractional n-forms can be constructed [|5(], pl| . For instance, the exterior 
derivative of the 1-form 

u = (dxTUix) (3.6) 
is d a co(x) = (d«f„)(x) (dx^) a A (dx u ) a . 

3.1 Interpretation of fractional gradients 

Fractional gradients admit a geometric interpretation [^, 60| , illustrated in figure ^| 



Already in one dimension, the fractional derivative of an a-differentiable function / 
on K can be expressed as 

(*«/)(*) = lim ( T °f)(y>*)-F»f)(x,x) (37) 

y-¥x+ (y - x) a 



where T a is a mapping suitably defined (see |5BL |5U] for details in the case of the 
Nishimoto derivative). In particular, 

(T a f)(x + h,x) = f(x) + d a f x (h) + h a e(h) , (3.8) 

where h > 0, d a f x (h) = h a (d a f)(x), and rim/ l _ s . £(^) = 0. (Related to this equation 



or variations on the same theme, one can develop a fractional Taylor expansion [82 



[j84}j.) In many dimensions, one can consider a directional fractional derivative and a 
bilinear mapping (T a f)(y, x) acting on MP <g) MP . Let (grad a /% := d^f be the /z-th 
component of the fractional gradient of a function /. Indicating a vector a; as a point 
M in R D , let the gradient exist at the point M G M D . If c?M is an infinitesimal 
vector displacement of the point, with coordinates (dx 1 , . . . ,dx D ), one can identify 
a fractional displacement dM a as the vector ((dx l ) a , . . . , (dx D ) a ). Therefore, the 
differential d a f is given by 

d a f = gmdj ■ dM a . (3.9) 

The set £ of points M satisfying (T a f)(M, M ) = c, where c is a constant, is called 
a level surface passing through Mq. The point Mo + dM in a neighborhood of Mo 
belongs to E, but the point M + <iM a does not. In fact, the latter determines 
another level surface (T Q /)(M, M ) = c', at an angle (3 with E determined by 
cos/? = |c — c'|/||(iM a ||, where || ■ || is the norm equipping the vector space. Since d a f 
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vanishes on £, the vectors grad a f and dM a are orthogonal, so that grad^/ is not 
orthogonal to (tangent vectors on) £ at M for < a < 1. The projection of grad^/ 
on the unit vector n normal to S at M has modulus grad Q / ■ n = |grad a /| sin/3. 
When a = 1, the two level surfaces coincide and c = c', (3 = 0. 





n 






grado/ \ 


M,i+clM a / 








/ - 1%+dM 







Figure 2: Geometric interpretation of fractional gradient. The symbols are explained in 
the text. 




3.2 Which calculus? 

We define fractional Euclidean space £® of order a as Euclidean space R D endowed 
with a set of rules Calc a = {d a , I a , . . . } of integro-differential calculus, a measure g a 
with a given support, a natural norm || • ||, and a Laplacian /C: 



(3.10) 



Different sets of fractional operators in Calc a can correspond to inequivalent frac- 
tional spaces. We should now make a commitment on the type of derivative and in- 
tegration operators acting in £®, and choose between left and/or right integrals, left 
and/or right Caputo or Riemann-Liouville derivatives, mixed operators, the range 
of a, and so on. In fact, there are many more fractional derivatives we could have 
listed here, so the choice is actually larger. For the purposes of this section, the only 
ingredient we need to specify is calculus. In particular, for finite- volume calculations 
we pick the left fractional integral with < a < 1, while we elect the left Caputo 
derivative d a with < a < 1 as the building block of the differential structure of 
£^ ■ We postpone the choice of the Laplacian /C to section |573| , after the discussion 
of the so-called harmonic structure of fractal sets. Furthermore, we anticipate that 
the measure of £® is, for each direction, the Weyl measure g a (x) = x a /T(a + 1), 
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with support on the positive real axis. Here we are interested in local properties of 
fractional Euclidean space, so we do not spell the reason why the support of g a can 



be chosen as [0, +00) and not something else. That is given in ||42 |. 
Let us justify in detail the calculus. 

• Fractional versus non-fractional calculus: Fractional calculus is among the 
most studied and best manageable frameworks generalizing ordinary calcu- 
lus and, to the best of our knowledge, there are almost no other examples of 
calculi mimicking fractal behaviours in physics and statistics such as anoma- 
lous scaling of dimensionality and discrete scale invariance [f|T|. Mathematical 



properties of fractals, however, can be analyzed by Connes quantized calculus 



851, which we shall not consider here. 



Left versus right operators: While left operators involve integration from an 
initial point Xq up to the arbitrary point x, right operators integrate over 
the complementary sub-interval In fractional mechanics, where x = 

t is time and one studies dissipative classical systems, the output of these 
operators should depend on the past (rather than future) evolution of the 
system, and left derivatives seem more natural. However, on one hand in 
our context we do not have this type of physical interpretation and, on the 
other hand, both types of operators appear in fractional Lagrangian systems 
because of eq. ( p.64| ). Fortunately, "left versus right" is a non-issue, since 
the two classes of operators are actually the same under coordinate reflection 
(eqs. ( |2.11| ), ( p. 12| ) and ( |2.23| )). Further discussion of this point can be found 



111 

< a < 1: This choice is empirical. One wishes to obtain a model of space 
whose dimension is smaller than the topological dimension D of the embedding 
space. For a natural isotropic distribution of fractional charge over the D 
coordinates, this is achieved precisely for this range of the order parameter. 
The range of a will be further restricted in section |3~3. 



Caputo versus Riemann-Liouville (and others): The Caputo derivative carries 
several advantages over the Riemann-Liouville operator. 

(i) First, the Caputo derivative of a constant is zero for any xq, while for the 
Riemann-Liouville derivative it is so only when xq = —00. This helps in 
rendering fractional calculus more akin to the ordinary one without affect- 
ing other properties. We have seen, for instance, that the fundamental 
theorems of calculus are generalized in a simple way only for the Caputo 
derivative. Another difference, fully appreciated only when doing tensor 
calculus on curved manifolds, is that the frame metric 77^ is indeed a con- 
stant matrix for the Caputo differentiation; hence, tensor calculus keeps 
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many of its usual rules. 4 In |42|], we shall see that it is possible to define a 
simple fractional generalization of Poincare symmetries precisely because 
of this property of the Caputo derivative. 

A popular reason to prefer Caputo over Riemann-Liouville operators is 
the existence, in the former case, of a standard Cauchy problem, where 
one needs to specify only the first n ordinary derivatives |45|| . On the 
other hand, the Riemann-Liouville operator requires the specification of n 
initial conditions of the form lim x ^, Xo (I k ~ a f)(x) , k — 1, . . . , n, which have 
no clear physical interpretation. Related considerations hold for general 
non-local theories, where the Cauchy problem must be reinterpreted; such 



a reinterpretation is known only for special non-local operators (e.g., [92 



(iii) There is, moreover, a simple but not very well known argument setting 
the Caputo and Riemann-Liouville derivatives on a different footing [j45|J , 
and favouring the former as the most plausible operator to appear in 
dynamical equations. Taking the a — > n limit of eq. (|2.15|) , one obtains 



n-l 



(d n f)(x) = Ud n f)(x) -Y^d n - l - j 5{x - x )(dff)(x ) . (3.11) 

This is nothing but the relation between the integer classical derivative, 
on the left-hand side, and the derivative in the sense of distributions on 
the right-hand side. Therefore, one can consider the Caputo derivative as 
the fractional generalization of classical differentiation, and the Riemann- 
Liouville derivative as the fractional generalization of functional differenti- 
ation. In this respect, the Caputo operator is a much more natural choice 
for the fractional derivative in actions defined on £®. 

(iv) Another difference between Caputo and Riemann-Liouville derivatives 



emerges within initialized fractional calculus |52|, P4}|~II98I1- We mention 
it for the sake of completeness, although it has no impact in the present 
discussion. Heuristically, while repeated use of ordinary integration and 
differentiation generates arbitrary integration constants which are fixed by 
the boundary conditions of the problem, fractional integration is a continu- 
ous operator leading to a non-trivial entire function, called complementary 
(or complimentary). This function is determined by additional input on 
the boundary conditions. More precisely, in classical fractional mechanics 
one is interested in observing a system starting from a time later than 
the initial time to when the system began to evolve; call f(t) the observed 



4 



This fact was recognized in studies of manifolds with non-holonomic structure 
framework employed to integrate non-linear dynamical equations such as Einstein's pC 
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history, a continuous function on [i*,t]. The effect of past history is in- 
corporated in the complementary function. It turns out that the inferred 
history for the Riemann-Liouville derivative is continuous throughout the 
evolution period, and it is described just by the function f(t), analytically 
continued to the whole interval [to,t]. On the other hand, if the Caputo 
derivative is required to be equal to the initialized Riemann-Liouville 
derivative, the analytic continuation of f(t) is the constant fit) = f(t*), 
for t < t* [p5[| ; therefore, derivatives of / are discontinuous at £*, and 
Caputo initialization effects seem not to be properly taken into account 
In our case this is not an issue, both because we do not ask the two 



derivatives to coincide (we simply make one choice and keep it all along 
the way) and because the context of fractional evolution is quite differ- 
ent. There is no connection between the lower terminal of integration, 
representing a boundary of space (time), with the notion of performing 
observations at space(time) points in the bulk. Here, all functionals are 
defined to have support in the domain of the fractional integrals, so there 
is no need to extend the physics outside this domain. 



v 



) The choice of measure (|2.85|) and (|3.1|) may seem too restrictive within 



all the possible choices of calculus, but this is not the case. On one hand, 
the more general arbitrary Lebesgue-Stieltjes measure ( |1 . 1|) is technically 
untractable [p3[ f40| , and not much progress can be done even for the abso- 
lutely continuous measure (|1.3| ) without specifying the profile v[x) |23|, |39|1 . 
On the other hand, there is no apparent advantage in choosing specific 
profiles different from ( 2.85Q and ( |3.1|) , which is highly anisotropic in the 



embedding coordinates. For instance, one might object that anisotropic 
measures do not have the symmetries we observe in Nature, such as ro- 
tation invariance, and one might advocate isotropic measures such as a 
power-law profile v(x) ~ |ai|~ D ( 1_Q ). Integro-differential operators gen- 



erating this type of measure constitute the so-called Riesz calculus |)46 
section 2.10], which is another version of fractional calculus. In principle, 
it would be possible to define fractional spaces also with this calculus. 
However, the factorized measure ( ft.l|) allows us to define a natural set 
of geometric coordinates (section |3.3| ) and a simple differential structure 
(section |3.4| ), while mixing coordinate components would result in some- 
thing more difficult to deal with. Furthermore, in fact, we do not need 
the measure to have any of the ordinary Poincare symmetries, because 
these turn out to be accidental symmetries of fractional spacetimes in the 



infrared fly, @2 



Different choices of calculus are not excluded. In principle, they will lead to differ- 
ent quantitative details (that is, in the mapping between embedding and geometric 
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coordinates and in the definition of the geometric integro-differential operators). Oc- 
casionally, we will compare the left theory with one with mixed operators. We saw 
that these operators do not share many of the simple properties enjoyed by each 
individual sector, and because of this we will not pursue a full description of a mixed 
theory. Nevertheless, an example of volume calculation will show how geometry 
changes from a single-sector to a mixed-sector formulation. 

It is important to stress that inequivalent fractional spaces are expected to share 
the same physical features. The reason is that the scaling property ( |3.2D of the 
fractional measure is unaffected by the choice of differential operators; see sections 



3.5 and 3.6 



3.3 Geometric notation 

The geometry of fractional space is characterized by a set of integro-differential oper- 
ators, specifying both the measure on the space and its boundary, and the notion of 
distance (shortest path) between two points. The latter can be and has been derived 
within fractional calculus, but these results are more transparent after introducing 
a novel notation, which we call "geometric" and also carries part of the physical 
interpretation of the model. 

The central idea is to regard a fractional manifold as embedded in an ordinary 
.D-dimensional manifold constituting a mathematical ambient space. A fractional 
manifold is defined by the multiplet of objects (|3.10|) but with more general embed- 
ding. By embedding we mean the ordinary manifold one would obtain by setting 
= 1 in the measure for all /i. A general embedding is a smooth manifold endowed 
with a metric g^. For instance, in the embedding defining fractional Minkowski 
spacetime M.® is ordinary Minkowski spacetime M D , with metric g^ u = r]^ u . In the 
case of fractional Euclidean space £®, the embedding is M D and g^ v = 5^ v . 

While /i = 1, . . . , D, are embedding coordinates, on the fractional manifold 
a natural coordinate system is provided by the "fractional" (or "geometric" ) coordi- 
nates g M , eq. ( |2.69| ). The symbol q leaves implicit both the coordinate dimensionality 
and the dependence on the embedding coordinate domain, but this does not differ 
substantially from what one does in ordinary geometry. In fact, hiding a fixed a 
in q is tantamount to hiding the scaling [x] = —1 in the symbol x. Similarly, one 
typically specifies a coordinate system {x} and its domain Dom(x) separately, and 
not as a joint symbol XDom(x)- The coordinates ( |2.69| ), sometimes called "general- 



ized", were introduced in the special case xq = (and without the T factor) in 



[ pq , |67| , |71| , |72| , jS0| , pl| in the context of dissipative mechanics. 

Equations ( |2.22|) and (|2.34|) state that there is no natural geometric coordinate 
for Liouville (and also Weyl) calculus, i.e., no function q(x) such that ood a q = 1. In 



fact, the right-hand side of eq. ( |2.69| ) diverges for xq = — oo (e.g., for global Carte- 



sian coordinates). This fact should not be of concern, since this is just a mapping 
stating how embedding and fractional coordinate systems are related to each other. 
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A singular mapping does not imply a pathology in the embedding coordinate system, 
and in fact the final expressions of geometric integrals are perfectly well defined in 
the language of fractional calculus, even for |sco,i| = 00 • I n other words, even if there 
are no geometric coordinates in pure Liouville and Weyl calculi, fractional spaces 
equipped with these operators are still meaningful. 

The "geometric" differential associated with the fractional coordinates q is just 

d a : 

d:=cP, [d] = 0. (3.12) 

We change notation to avoid confusion between the label a and space indices //, u, 

Prom eq. (|2.70|) , 

d = dqd%, dq = {dx) a , (3.13) 

where the "geometric" derivative is the Caputo fractional derivative regarded as a 
function of q, 

d« := A = d: - (3.14) 



in particular, d®q = 1. Note that eq. (|3.13|) can be taken as the definition of q 
via d a q also in Liouville calculus (xq = —00; compare eq. (2.9) in |99]|). Also, 
dg 7^ d q = (dq/dx)~ 1 d x : comparing with eq. ( |2.62| ), 



a 



df 1 

+ -: 



f(x) ~ f{x ) 



r(a)r(2- a) dq gr(l-a)r(l + a) 

sin[7r(j - a)} T(l + a) , ,, a(r)iru 



i=2 



n(j-a) T(l+j) 



(3.15) 



In many dimensions, 



d:=ikf <9" 



0° 



d a q^ 



(3.16) 



By virtue of the unique property d®l = typical of the Caputo fractional derivative, 
d^q u = 8" . The symbol will indicate both the partial fractional derivative with 
respect to 2^ and the one with respect to g 1 '; the context should make the distinction 
clear. 

Finally, we define the "geometric" integral 



T(a) 



.Co 



dx' 



x — x 



l-a 



0. 



(3.17) 



The symbol in the left-hand side of the first equation is borrowed from the standard 
notation for mean integrals, which are, of course, out of the present context. The 
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right-hand side was first introduced by Tarasov [50]. A definite integral over the 
whole interval is simply obtained by setting x = x\. In one embedding dimension, 



dq' 



(dx 



/\l—a 



T(a) J XQ (.'• - r 



l\l-a 



{dx') 



r{i} 

(X-X ) c 

T(l + a) 



Therefore, posed f(g) = f(x), eq. ( j2.60| ) is equivalent to 

f di(q') = f(q) - f(0) , 
Jo 



(3.18) 



stating that geometric integration in fractional coordinates has formally the same 
properties as the ordinary integral. By formally, we mean that there will be a differ- 
ence in the functional space over which the integral operators act. 

Geometric integrals with different lower extrema are obtained either by a initial- 



ization prescription |94[ or by the composition law 

- x o y 



q*,i = 



r(i + a) 



< q* < qi 



' q* JO JO 

Multiple integrals in the coordinate system /j, = 1, . . . D} follow through, 



jd D q := j dq 1 ■ ■ ■ j dq D . 



(3.19) 



Each fractional coordinate g M is mapped into an embedding coordinate x^ with a 
given domain. Examples of nested integral will be seen in section |3.5| . 
We summarize the three notation systems employed so far in table |]. 





Geometric 


Fractional calculus 


Stieltjes measure 




formalism 


formalism 


formalism 


Coordinates 


q 


(x-x ) a /T(l + a) 


X 


Measure 


dq 


{dx) a 


dg(x) 


Integration 


f 


r 


I 


Differentiation 






d 



Table 1: Equivalent formalisms describing calculus on a fractional manifold. The 
Lebesgue-Stieltjes measure formalism, with generic measure g, is the most general but 
it is often impractical to perform calculations. 
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3.4 Metric and distance 

In the language of first-order general relativity, arbitrary coordinate transformations 
define frames which, in turn, determine the metric. Therefore, the notion of line ele- 
ment naturally emerges. First-order formalism is somewhat an overkill when dealing 
with Euclidean space, but one can foresee obvious applications in more general sce- 
narios. 

Consider two coordinate systems {x 1 } and {y^ 1 }, the first (denoted with capital 
Roman indices) being the Cartesian system and the second a generic curvilinear one. 
The exterior derivative of order a can be written in both systems as 

{dx I ) a df = d a = (dy») a d« . 

Applying eq. (|2.68|), we get 



(dx J ) a = {dyT ^M =■ (dyTet , (3-20) 

r (l + of) 



where erj, a D x D matrix, is the fractional generalization of the vielbein. To avoid 
confusion with space indices, we shall omit labels a for fractional vielbein and metric, 
using the same symbols as in ordinary space. 

Some example of coordinate transformations are given in for Riemann- 
Liouville fractional calculus. We can easily give the general form of the fractional 
Jacobian J a in a calculus of order a. Let J{y) = \dx(y)/dy\ be the Jacobian of 
a coordinate transformation from the system of coordinates x to y. The fractional 
Jacobian J a is simply J times the ratio of measure weight factors: 

a- I 1 D\ v a{y l ) V a (y D ) 1 D /qqi\ 

*'<»•■■■'* )= ^m^'''^m J{v '-' v >■ (3 ' 21) 

For instance, for left integrals with upper terminal 

v a (y l ) v a (y D ) 



v a [x l {y)] v a [x D (y)] 



11 DDI 1~ Q 

yi -y y* - y 



xl(y*) - % l (y) x?(y*) - % D (y) 



In compact notation, the integral measure transforms as 

dg a (x) = dg a (y) J a (y) , J a (y) = V ™ V J{y). (3.22) 

As in ordinary calculus, the measure in the new coordinates does not factorize over 
the directions due to the non-trivial coordinate dependence of the fractional Jacobian. 
Expressing the 1-form ( ^.6| ) in y coordinates and reversing the transformation, 

lo = (dx J rfj(x) = (dyTe J Jj[x(y)} = {dx'f ^e J J \j[x 
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from which it follows the relation 



eH = 5j . (3.23) 



One can also define the fractional metric 

9nu := Vueiei , (3.24) 

where rju = Sjj is the Kronecker delta in Euclidean space. In turn, the fractional 
metric gives the fractional line element 

ds a := [g^{dx^) a ® {dx u ) a Y* , (3.25) 

or, in geometric notation, 

ds 2 = g^Aif ® dq v . (3.26) 

This result for fractional two-forms suggests a natural definition of the distance 
between two points. The metric of fractional Euclidean space £® is g^ u = S^, so in 
geometric coordinates it is expressed as 




A a (q, q') := y/A^A^ = [A a (q\ g' 1 )] 2 + ■ ■ ■ + [A a (q D , q'°)f . (3.27) 

This only resembles an ordinary Euclidean distance, since A a (q, q') ^ \q — q'\. From 
eq. ( |3.25| ), the coordinate distance is the 2a- norm 

(D \ ^ 

E W - y"\ 2a J , (3.28) 

where A(x Al ,?/ M ) = |x M — This is a norm only if a > 1/2, i.e., when the triangle 
inequality holds. Therefore, we can further restrict a to lie in the range 

(3.29) 

From the perspective of differential forms, the 2a-norm is the natural distance in 
fractional space, which is a metric space. In fact, one should not confuse eq. ( |3.28j ) 
with the choice of a p-norm (all topologically equivalent) in a given space: as a 
changes, so does the geometry of space. 

In a generic fractional geometry with a / 1, the Pythagorean theorem is not 
valid and the shortest path between two points is neither a straight line (Euclidean 
distance) nor unique. The case a = 1/2 corresponds to the so-called "taxicab" or 
"Manhattan" distance, given by the rectilinear distance along the axes. In D — 2, 
circles in this geometry are diamonds with edges at 45°; the inclination of the edges 
is fixed, taxicab distance not being rotation invariant (figure ||). As a increases from 
1/2 to 1, the faces of the diamond become convex until they merge into an ordinary 
circle. These figures are called supercircles (a particular case of superellipse, or Lame 
curve). See figure f|. In D = 3, taxicab spheres are octahedra (figure |5|). 
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Figure 3: 1-norm and taxicab geometry in two dimensions, (a) Left panel: the shortest 
path between two points is not unique, (b) Right panel: circles of radius R are diamonds 
with edges at 45° with respect to the coordinate axes, \x\ + \y\ = R. 




-1.0 -0.5 0.0 0.5 1.0 0.0 0.5 1.0 1.5 2.0 

X q\ 



Figure 4: Left: circles in two dimensions with unit radius in 2a-norm, \x\ + \y\ = 
1; increasing thickness corresponds to a = 1/2, 3/4, 1. Right: the same circles in left 
geometric coordinates q± = (x — xo) a /T(a + 1) and q2 = {y — yo) a /T(a + 1), with xq = 
-1 = Vo- 

3.5 Volume 

We have already described fractional operators in one dimension, and their repli- 
cation to D copies is straightforward. As one might expect, there is a fractional 
analogue of areas and volumes, but their scaling properties and values will differ 
from their ordinary Euclidean counterparts. 

Let M. C be an arbitrary region in space. The fractional volume of A4 is 
(e.g., H) 

D 

where each fractional integral is denned on an interval [ These intervals 

can be always chosen so that they cover Ai. If Ai is a rectangular region, then all 
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Figure 5: Spheres in taxicab geometry are octahedra, \x\ + \y\ + \z\ = R. 

the coordinates Xq and x^ are constant. Otherwise, one can compute the volume Vm 
as a nested integration over elementary domains, as in ordinary calculus. 

We give an example of multiple integral in the calculation, in D = 2, of the 
volume V^ 2 \R) of a 2-ball with radius R and centered at the origin. 5 With the 
a-norm distance (|3.28|), this is the area enclosed by a supercircle B2, i.e., the set 



B 2 = {(x,y) : \x\ 2a + \y\ 2a <R 2a } . (3.31) 

When a — 1, this is a disk with radius R; when a = 1/2, it is a diamond with 
vertices (±-R, 0) and (0, ±R). We integrate first y from —r(x) to r(x), where 

r(x) = (R 2a - \x\ 2a ) 1/(2a) . (3.32) 

Then, we integrate in x from —R to R: 

V {2) (R) = I^ R {^ r{x) , r{x) {l}} 

] [2r(x)]° \ 

R 2a , (3.33) 



where 



/2,«= p^T f dx {I - XT' 1 

1 (a) J-i 



2a 



1 /"I , 

= / dx Vl-x 2a [(l - x) " 1 + (1 + x)"- 1 ] . (3.34) 

We will reconsider this prefactor later for general D. Other examples of double 
integrals can be found in |50fl . Working in geometric notation, the two fractional 



^Clearly, the final result will not depend on the location of the disk in the coordinate plane. 
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coordinates are 



a = £±J9! 6P , R] , W ==^L, (3.35) 

1(1 + 0) 1 (1 + a) 

, 2= J|+rM! eM , r-^M. (3.36) 

1 (1 + a) 1 (1 + a) 



Therefore, 

y( 2 )(ft) = / dgx / dg 2 = / dfc f( gi ) = 2^T(1 + a)/ 2 , a # 
Jo Jo Jo 



2 

) 



coinciding with eq. (|3.33|) . 



As we said, different presentations of the fractional operators lead to inequivalent 
fractional spaces, where areas and volumes are weighed according to the measure in 
the integral. For example, the measure weight of the a = 1/2 unit disk is v(x,y) = 
[(1 — x)(l — \x\ — y)] -1 / 2 , which is heavier for points (x,y) ~ (1,2/) and (x,y) ~ 
(x, 1 — |x|), corresponding to the right vertex and to the upper edges of the diamond. 
This is depicted in the left panel of figure |6|. Using right integration J a [x]J Q [?/], the 
measure is v(x,y) = [(1 + x)(l — \x\ + y)]~ l l 2 . However, due to the symmetry of 
the object (eq. ( |3.34j ) is invariant under x — > —x), the value of the area is the same; 
this is not true in general. On the other hand, for mixed integration J a [x]/ a [j/] with 
(£*,2/*) = (0,0), we have 

= I-R,0,R {2Io,r(x){l}} 

2[r(x)] c 



TO 

1 -R,0,R 



T(l + a) 



2<i 



r(i + a 

where we used the fact that the integrand is even and 



R 2a , (3.37) 



ha = [ dxx^VT^x^. (3.38) 

r(«) Jo 



Therefore, in this case the measure weight for a = 1/2 is v(x,y) = |xj/| _1//2 , which 
diverges along the axes x = and y = 0; see the right panel of figure || 

The angle factors / Qj2 and / 2 ,a represent the ratio of the area of the disk and a 
certain power of its radius. Measurements of both would then determine the angle 
factor and, in principle, provide an experimental discrimination between left/right 
and mixed theories. Let us recall the overall physical picture outlined in section |L2l 



- 40 - 




Figure 6: Density plots of the area measure of a 1-norm disk. The integration measure 
weight is represented in light to dark shade, darkest shade being points where it diverges. 
From left to right: left integration, right integration, mixed integration. 

The fractional space constructed in this paper has fixed dimensionality and there is no 
dimensional flow. So, the dimension is always different (possibly very different) from 
D at all scales. To get some physics, one must generalize to spaces which are multi- 
fractional and whose dimension changes with the probed scale. This is overviewed in 



[] and done in [42]. In the full theory of [42], spacetime has dimension close to D at 
large scales, and local geometry can be tested against deviations from the ordinary 
one. With this setting in mind, one can envisage local measurements of geometry in 
vacuum, made at scales and in conditions where gravitational effects (such as tidal 
forces) are negligible, to check if geometry in a local inertial frame is Euclidean. 
For instance, one could take experiments on the equivalence principle and Lorentz 
violation and, by reverse engineering, place experimental bounds on the angle factor. 
In standard theories, at human scales (say, below 1 km) Euclidean geometry holds 
and the theoretical value of the angle factor is ir. However, if a fractional theory with 
dimensional flow was a correct description of Nature, we would expect the parameter 
a not to be exactly equal to 1 at sufficiently large scales. In other words, an expansion 
in a = 1 — e/D, 0<e^C 1, would yield a correction to the Euclidean angle factor, 
which can be constrained from above by experiments |32[ . 

We can find the form of this correction in arbitrary dimension. To begin, we 
prove by induction that the volume of a closed D-ball 



B 



D 



D 



■ £ 



x 



fi 1 2a 



< R 



2a 



(3.39) 



is 



V {D) (R) = Q Dia R 



Da 



(3.40) 



where £ln,a is the volume of a unit ball. We first work with fractional calculus and 
then in geometric notation, to show how the latter is more transparent. Suppose 
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eq. (|3.40|) true in D — 1 dimensions for a ball of radius r(x) = (R 2a 
Then, integrating also in the direction x, 



1 2.120^1/(20) _ 



V {D \R) 



x=x'/R 



n 
n 



D-l,a*-R 



')} 



D-l,a 



Da 



(3.41) 



where 



1 



D,a 



r(a) 
1 



dx (1 — x 



,a-l, 



1 — \x 



2a\ 



-1 
1 



da; (1 - x 2a ) — [(l - x)*- 1 + (1 + x 



(3.42) 



In geometric notation, the fractional coordinate over which one integrates is given 
by eq. ( p.35|) . In terms of the fractional radius 71, the the volume (|3.40|) scales as 



I'D, a 



(3.43) 



Then, given a (D — l)-ball of radius r = (2r) Q /r(l + a), one has 



v {D \n) 



n 



dqV {D - l) [?(q)} 
u D -i, a J <dqr D 1 (q) 



(3.44) 



= 2- a T(l + a)u D ^J Dta n D , 

in agreement with eq. (|3.41|) . Formally, this is the same calculation as for (one orthant 
of) a .D-ball in Euclidean space, the difference being in the angular factor. Thus, 
from the point of view of the observer in the fractional manifold, the ball scales as it 
should. However, fractional coordinates/distances have anomalous scaling (x — > Xx 
implies q — > X a q), so the embedding volume of the ball scales as R Da . When D = 2, 
one recovers the explicit result for the superellipse. 

The angular factor Qd,o is obtained by solving the recursive equation flo, a = 
&D-i,afD,a with initial condition f2 lQ , = f la = 2 a /r(l + a): 

D 

£l D , a = l[f nia . (3.45) 



n=l 



We were unable to compute Qd,o explicitly in finite form for general a, but one 
can do so for several special cases. For instance, when a = 1/2, one obtains 
eq. 3.197.3] 



When a — 1 — e/D, < e < 1, we get the standard result plus corrections. In fact, 



DA- 



e 

D 



dx (1 — x ) 2 
l 

axil — x ) 2 



o 



2 7 + ln(l - x 2 ) - (D - 1 



x 2 hax 
1-x 2 



21 



0(e 2 



where 7 = —"0(1) ~ 0.577 is Euler constant and i/j is the digamma function. Using 
formulae 3.251.1 and 4.253.1 of 



^r(g±i) r e 



27 + 



D + l 



1>[l 



Then, 



0(e 2 ) . (3.47) 



D 



&D,1- 



II ./;,, 



n=l 



27-V^| 3) + 



71=1 v 



n + 1 



where 



D 



n^i = n 



7T 



D/2 



0(e 2 ) , (3.48) 



(3.49) 



J r(f + i) r(f + i) 

is the standard unit volume. 

In tables and |] we compare these unit volumes with the standard ones, and with 
the result which would come from a traditional dimensional regularization procedure, 
where the topological dimension is formally expanded as D — e: 



£lo-e,i — ^d,i { 1 — - 



111 TT - (.'(— + 1 



+ 0(e 2 ). 



(3.50) 



D 


£1d,1/2 


£1d,1/2 


2 




v/2 


^7T = 4 

7T 


3 


(V2- 


» 0.88 


8 4ir 
7T3/2 3 ~ ° 


4 


2 v / 2(v / 2- 


- 1) (1 - % ) » 0.43 


16 TT 2 
IT* 2 



Table 2: Volume of unit D-balls in various dimensions for a = 1/2, in left/right and 
mixed theories (with x* = 0). 

Left fractional expressions are not as neatly symmetric as those of integer calculus 
because fractional integral operators are not even in the embedding coordinates. In 
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the mixed theory, the volume of a D-ball scales as in eq. ( [3.40 ), but with angle factor 
given by 



D 



^D,a — Y\ fn,< 



n=l 



where 



r a) 



dx x a (1 — x 



2a 



D-l 
2 



(3.51) 



(3.52) 



and fli ta = f\ t a = 2/r(l + a). This integral can be expressed in terms of T functions 
and one has 

QD > a = [T(l + a)} D • (3 - 53) 

In particular, for a = 1/2 the volume is enhanced by a factor (2/y / 7r) D , and in 
general it is considerably greater than in the left theory (table ||). Expanding in 
a = 1 — e/D, one has 



S7 D ,i < =fi A1 [l + e(l- 7 )] + 0(e 2 ). 



(3.54) 



The coefficients in the e corrections have opposite sign with respect to the left/right 
and dimensional-regular izat ion cases (table |3|), and are one and the same for any D. 
As a consequence of eq. ( |3.53| ), which define the effective constant tt := ir[l + e(l — 
7) jn\ , these corrections can be written relative to it and are the same for D = 2n and 
D = 2n + 1 dimensions in the mixed theory. Anyway, e is expected to be constrained 
to extremely small values by experiments, so the actual coefficients in front of it 
are not important except for highlighting this concrete comparison of inequivalent 
theories. 



D 








2 
3 
4 


tt(1 - 0.42e) 
-0.54e) 
^(1 -0.63e) 


7r(l + 0.42e) 
f (l + 0.42e) 
4(l + 0.42e) 


tt(1 - 0.36e) 
f (1 -0.22e) 
^(l-0.11e) 



Table 3: Volume of unit D-balls in various dimensions, for a ~ 1, in left/right and mixed 
theories (with =0). The corrections in traditional dimensional regularization are shown 
in the last column. 



3.6 Hausdorff dimension of space 

When dealing with exotic sets, it is important to define a sensible notion of dimension. 
Sometimes, the imprecise name "fractal dimension" is used to indicate one or more 
among the many possible (and inequivalent) definitions of dimension, which may 
create much confusion. Here we specialize to one such definition, the Hausdorff 
dimension. Before doing so, we make general remarks on dimension counting [|100| . 
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Let T be an object living in a D-dimensional space. To measure its volume, one 
can take the minimum number N(5) of n-balls, n < D with radius 5 centered at 
points in J 7 and such that they cover T (i.e., each point in J 7 lie in at least one ball). 
The number N(S) increases as 5 decreases, approaching the behaviour N(8) ~ 5~ dB 
as 5 — > 0. Then, the number 

rf B :=-hm^P (3.55) 
s^o In 5 

is called the box-counting dimension of J 7 . (Strictly speaking, there exist a lower 
and an upper box-counting dimension, given by the liminf and limsup, respectively; 
when they coincide, they reduce to eq. fl3.55| ).) For instance, if J 7 is a square or a 
disk, the covering of 2-balls will show that d^ = 2, if it is a cube or a 3-ball, a 3-ball 
covering will give d-Q = 3, and so on. Very irregular or fractal sets will not be smooth 
and their dimension will be, in general, non-integer (although there exist also fractals 
with integer dimension). Intuitively, a set with many irregularities will require more 
balls for being covered, and their number will increase faster than expected; a typical 



example is an irregular porous surface, for which d-& > 2 ||101|| . On the other hand, a 
surface with "too many holes" may require less balls than a smooth one. For a given 
n, the shape of the covering sets is not important, and one could use, for instance, 
n-cubes of edge length 5 instead of n-balls; what matters is how the volume of the 
probe scales with its size. 

The box-counting dimension above is only a particular definition of dimension, 
and it often proves to have a number of inconvenient properties. It is desirable to 
have a different (not just more general) notion of dimension as follows. The idea is 
similar, namely, to define a "minimal" covering for J 7 C 1R D , but now taking; covering; 
sets of different size. Let \U\ = sup{A(a;, y) : x, y e U} be the diameter of a set 
U C ~R D , i.e., the greatest distance A(x,y) between two points in U. A 5-cover of 
J 7 is a countable or finite collection of sets {U} of diameter at most S that cover J 7 : 
J- C U i Ui, with < |Z7i| < 5 for all i. If s > is a real non- negative parameter, one 
can define 



{2> 



^(jr) := ii min f <> y \Ui\ s : {£/*} is a 5-cover of 7 ) . (3.56) 

This limit exists (it can be also and +oo) and is a measure, the s-dimensional 
Hausdorff measure of J 7 . One can check that it is proportional to the n-dimensional 
Lebesgue measure for integer s = n < D (length, area, volume, and so on). The 
Hausdorff measure obeys the scaling property (in coordinate notation) 

(&(Xx) = \"l&(x), (3.57) 

where A > is the scale factor of a dilation x — > Xx. One can show that is 
non-increasing with s and there exists a critical value of s at which the measure 
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jumps from +00 to 0. This is the Hausdorff dimension (or Hausdorff-Besicovitch 
dimension) of T [102]: 



d n {F) := inf{s : g s {F) = 0} = sup{s : g s {F) = +00} . 



(3.58) 



This definition allows one to calculate o?h via the behaviour of the Hausdorff measure. 
The latter diverges for s < g?h, is zero for s > dn, and may be 0, +00 or finite at 
s = du- In general, 4 ^ 4^ 

Taking balls as the covering sets Ui, one defines a measure which jumps at the 
same critical value g?h of the Hausdorff measure ||100| . This determines a local, 
operational definition of the Hausdorff dimension g?h of a smooth set such as an 
integer or a fractional manifold of topological dimension D: dn is given by the 
scaling law for the volume of a D-b&ll of radius R: 



V {D \R) ocR d " 



Therefore, 



dn 



lim 

5^0 



lnV^) 



In 5 In 5 

where in the last step we formally expressed the volume as the measure of a D-ball. 

Thus, in section |3~5| we implicitly proved that the Hausdorff dimension of isotropic 
fractional Euclidean space is 



lng[B D (5)] 



(3.59) 



(3.60) 




(3.61) 



Also, we calculated volume corrections for a nearly integer dimension, when a ~ 1 
and 

d H = D-e. (3.62) 

The scaling of V^ D ' reproduces the estimates of fractal distributions in the fractional 
continuum approximation |66|, [J7], [H], [72], |103|| - ||T06| as well as the heuristic scaling 
of general Lebesgue-Stieltjes measures in certain regimes [^3], [39], |40[| . 



Since fractional space is smooth, one expects most of the inequivalent definitions 
of fractal dimension to collapse one into the other. For instance, the Hausdorff 
dimension of the product T\ x Ti of two fractals is greater than or equal to the 
sum of the Hausdorff dimensions of the two sets, but it is strictly equal if the upper 
box-counting dimension coincides with the Hausdorff dimension for either J-\ or J-" 2 
[ p.0U| , corollary 7.4]. The box-counting dimension of a fractional ball is the same as 
its Hausdorff dimension, 

d~B = dn , (3.63) 



hence the result ( |3.61| ). In section |4.2| , we will reobtain eq. ( |3.61| ) two more times by 
symmetry arguments of fractal geometry. 
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The class of fractional spaces is bounded by two limiting cases. When a = 1 
in all directions, one recovers D- dimensional Euclidean space, where the measure 
weight is uniform and all points are on an equal footing. A particle has full memory 
of its past history and the dynamics is determined by certain initial conditions. On 
the opposite side, when a = the measure weight is peaked at the boundary of 
space, but integration reduces to the identity operator. The action is defined at a 
point and the space is zero- dimensional. This "Pointland" universe has no memory 
whatsoever of its history and, because it has no extension, it has no dynamics at all. 



4. Fractional versus fractal 



As soon as fractals made their appearance in the literature |[L07|| , the anomalous 
scaling of fractional measures induced the perception that certain phenomena with 



fractal properties might be described by fractional calculus ||108| . Later, it was argued 
that the fractional charge a is related with the Hausdorff dimension of certain fractal 
objects |S3|, |109| |. Fractional equations can approximate, in some sense to be made 
precise, self-similar deterministic (also known as nested) fractals such as the Cantor 
set ||110 l and von Koch curves ||109|| , and random fractals such as (the trail and graph 
of) Brownian and fractional Brownian motion [ 109|| . Criticism on the results of ||110| 
and on the connections between fractional calculus and fractals [II 1 II, 11121 led to 



their clarification for the Cantor set ||113| , |114|| and their progressive generalization to 
self-similar sets (finite or infinite) generated by linear mappings (random self-similar 

8], |113| , |114j| , and generalized self-similar sets 



[[113| , |114|| , self-similar deterministic 
WW), generalized cookie-cutter sets where only the first similarity is linear [|1 1 5|, |116 



generalized net fractals (defined by contractions) where only the first mapping is 



linear ||117|| and, finally, generalized net fractals generated by non-linear mappings 

EMEU- 

By now, it is established that fractional systems are not indiscriminately equiva- 
lent to fractal systems. On one hand, there are features of deterministic fractals which 
are not reproduced by the simplest fractional systems. On the other hand, random 
fractals are indeed describable by fractional tools (section |4.4|) . Dynamical systems 
with fractal properties in certain static regimes have been successfully modeled by 
chaotic and fractional systems jK], 52[ . Fractional differential equations (such as 
the generalization of the Fokker-Planck-Kolmogorov equation) well describe, for in- 
stance, self-similar dynamics, Levy flights, and anomalous diffusion in chaotic Hamil- 
tonian systems and systems close to thermal equilibrium []64|, [66|, [H], [72], |121|| - fl"3~2 
(see ||133| , |134|| for reviews) . Fractional Brownian motion ||100|| is also related to frac- 



tional calculus. Fractal domains characterized by a mass distribution or correlation 
functions with anomalous scaling are, in general, very irregular at small scales; this 
is the case in many physical systems such as porous materials, colloidal aggregates 
and branched polymers ||101|, |135||-|]THH[. These media, however, can be considered 
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as continuous at scales much larger than the characteristic size of the irregularities, 



such the pores in porous media [ 103 1 . Fractional systems, therefore, can be regarded 



as continuum approximations where the detailed microscopic structure of these ma- 



terials is smoothened without loosing anomalous scaling [ 103 , 105 , 106 1. This picture 
holds not only for mass distributions but, for instance, also in the description of the 
propagation of electromagnetic waves in dielectric media | |139| , and in other applica- 
tions (e.g., ||106|| and references therein). 



To understand in what sense fractional models describe fractals, and whether 
fractional Euclidean space £® is a fractal, we first fix the rules of the game, analyze 
what properties a fractal should have, and compare these properties with those of £®. 
This is a natural starting point where to draw a more precise comparison. Perhaps 



the most universal qualities of fractals are [|10C 

1. a fine structure; 

2. an irregular structure; 

3. self-similarity. 

Despite the fact that there are counterexamples of fractals not possessing one or 
more of these features, one must rely on them rather than on a sharp mathematical 
definition. As a matter of fact, there does not exist a unique definition of "fractal," 
other than "I know one when I see one" 111401 . 



4.1 Fine and irregular structure 

A fractal J 7 has a fine structure if it has detail at every scale. Intuitively this means, 
first of all, that one can zoom indefinitely into a fractal and always meet points 
belonging to J 7 , and, secondly, that in doing so one will always see non-trivial details. 
A smooth manifold M. as well as fractional space £® can be zoomed in indefinitely 
(they are continuum structures) but they lack details at all scales. 

Typically, fractals are also too irregular to be described with traditional geo- 
metric tools. This means that ordinary calculus does not apply to very discontin- 
uous sets, and one must resort to rather advanced techniques to define measures, 
Laplacians, spectral theory, and so on. A smooth manifold M. does not satisfy this 
property. Fractional space does by definition, although this is associated with 
an asymmetry of measure weights rather than manifest irregularity. At this point 
it becomes clear why fractional models are regarded as approximations of certain 
fractals: they do possess properties 1. and 2., but in a rather "dull" way. So £® is 
a fractal, technically, albeit of a rather uninteresting type as far as these properties 
are concerned. 

This is true only for fractional Euclidean space, where a is real and fixed. When 
lifting both these assumptions, and allowing a to be complex-valued and vary with 
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the scale, it turns out that the structure becomes extremely rich, and much closer to 
that of genuine multi-fractal sets [42 . 



4.2 Self-similarity and self-affinity 

Many fractals are self-similar, either exactly, approximately, or statistically. Roughly 
speaking, a set J 7 is exactly self-similar if it is made of iV(A) copies of itself of scale 
A. Then, the similarity dimension or capacity of the set is ||141] , |142| 



For instance, a hypercube in D dimensions can be thought of as the union of N = 2 D 
copies of itself, each of size A = 1/2 with respect to the original. Then, dc = D. 
One can take finer subdivisions in N(X) = (1/A) D copies of size A, and obtain the 
same result. In this sense, ordinary Euclidean space is trivially self-similar: A can 
be chosen arbitrarily. Genuine fractals have a more interesting self-similarity: for 
instance, the von Koch curve is the composition of 4 copies, each 1/3 of the original, 
hence dc = In 4/ In 3; the Cantor set is made of 2 copies of the original, each scaled 
A = 1/3, hence dc = In 2/ In 3; and so on. In most cases, like those just mentioned, 



the similarity dimension coincides with the box-counting dimension (|3.55|) and the 



Hausdorff dimension ||100| , section 9.2]. It is instructive to see this via the rigorous 
definition of self-similar sets. The ensuing calculations, employing the technique 
of contractions so often used in fractal analysis, will be more lengthy than those 
stemming from ( |4.1|) , but they will also allow us to obtain very precise informations 
about the structure of fractional space and its Hausdorff dimension. 

We begin with sets in M. D . Consider a set of N maps Si : M, D — > M D , i = 
1, . . . ,N >2, such that 

A[Si(x) , Si(y)] < AjA(x, y) , x,y eR D , < A, < 1 , (4.2) 

where the distance A between two points is A(x,y) = \x — y\ in ordinary integer 
geometry. Any such map is called contraction and the number A; is its ratio. If 
equality holds, Si is a contracting similarity or simply a similarity; if, moreover, 
Xi = 1, it is an isometry. Thus, a similarity transforms a subset of R D into another 
set with similar geometry. Many fractals are invariant under contraction maps and 
can be expressed as the union 

N 

F = {JS t (F). (4-3) 
i=i 

Given N contraction maps Si, T exists, is unique, non-empty and compact. The 



writing (O) is an "embedding" definition of a fractal, which is thought of as a 
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subset of Euclidean space. There is also a topological presentation which requires no 
embedding. For any non-empty compact set U one can define the transformation 



N 



S(U) :=\JSi(U) 



(4.4) 



i=l 



and its fc-th iterate S k := S o ■ ■ ■ o S. If Si(U) C U D J 7 for all z, then one can show 
that 



.F=ns*(co 



(4.5) 



fe=i 



In practice, this means that a fractal can be constructed by iterations of contractions. 
The k-th iteration may be regarded as a pre-fractal, an approximation of J 7 . The 



Hausdorff dimension of these fractals is bounded from above |100 

TV 

dn{F) < s , where Af = 1 . 



(4.6) 



When the <Sj in eq. ( |4.3| ) are similarities, the attractor J 7 is called an exactly (or 
strictly) self-similar set, and it is a union of smaller copies of itself In general, 

one further requires that the similarities obey the open set condition: namely, there 
exists a non-empty open bounded set U D T containing a disjoint union of its copies, 



N 



UD{JSi(U). 



(4.7) 



i=i 



If this condition holds, then one can prove that 




(4- 



A heuristic quick proof of this formula makes use of the scaling property ( j3.57|) of 
the Hausdorff dimension. For a self-similar set (4.3), 



N 



TV 



^(^) = ^^[^(^)] = ^A^(^) 



(4.9) 



t=i 



i=i 



and assuming that ^(J 7 ) is finite at the critical value s = c?h(^ 7 ), one can divide by 
^(J 7 ) to obtain flOP . 

Thus, if one can define a set J 7 via similarities, eq. ( |4.8|) gives the Hausdorff 
dimension of J 7 . Let us reconsider some of the examples mentioned at the beginning 
of the section. The middle-third Cantor set is defined by 



3 X ' 



<S2(x) = oX + 



2 . 

3 ' 



(4.10) 
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then, 1 = 2(l/3) dH implies g?h = In2/ln3. Another self-similar fractal is the 
Sierpihski triangle or gasket, the attractor of three similarities of ratio 1/2; there, 
1 = 3(l/2) dii implies g?h = In 3/ In 2. Other fractals, such as the non-linear Cantor 
set, can be defined by non-linear contractions or similarities. A more trivial but, for 
our purpose, instructive case is T = 1R D . In fact, one can work instead with a com- 
pact subset, the unit hypercube T = [0, 1] D = [0, 1] x • • • x [0, 1], where the extrema 
are chosen without loss of generality; then, the natural unbounded extension of the 
set will have the same Hausdorff dimension. Take first the case D = 1 and the interval 
T = [0, 1]. This can be expressed as the union [0, 1] = [0, A]U[A, 1], where < A < l. 6 
But this is equivalent to define two similarities such that J 7 = S\(J^) U S2 {J 7 ), where 

Si(x) = \x , Ai = A , (4.11a) 

S 2 (x) = (1 - X)x + A , A 2 = l-A. (4.11b) 

These maps satisfy the open set condition, as one can verify simply by taking open 
intervals. Therefore, 

1 = A dH + (1 - A) dH d H = 1 ■ 

In more than one dimension, it is convenient to set A = 1/2. In two dimensions, one 
has a square given by four smaller copies, and 4(l/2) dn = 1 yields g?h = 2. In general, 
a hypercube can be expressed by N = 2 D similarities (as many as the number of 
orthants) with equal ratios Aj = 1/2, so that 

1 = 2D (£) H ^ d n = D. (4.12) 

~R D is not a genuine self-similar fractal because it is not defined by similarities with 
a fixed < A < 1: the similarity ratio A is arbitrary. Also, a finite number of 
iterations S k (actually, just one, k = 1) is sufficient to obtain the set, which coincides 
with its pre-fractal approximation. In this last sense, Euclidean space is not a non- 
trivial self-similar fractal. Furthermore, M. D is not completely characterized by the 
similarities (|4.11|) . In the above construction for a unit hypercube, the similarity 
iSi(x^) = Ax M + a is just a contraction and translation equal along all directions, 
but M. D enjoys many more symmetries, including inequivalent contractions/dilations 
along different directions, rotations (A[S(x),Xo) = A(x,xq)), translations (S(x^) = 
x^ + a M , where also a now is a vector) and reflections (Six 11 ) = a M — x M ). These are 
examples of affine transformations, linear mappings of the form 

x" 1 = S(x") = A»x u + a", (4.13) 

6 The argument can be carried out verbatim for the semi-open interval [0,1), where Si(T) n 
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where A is a D x D matrix. Linear similarities are a particular case of affine trans- 
formations. The attractor of a sequence of affine transformations is a self-affine 
set: 

N 

F = {JlMF) + ai ]. (4.14) 
i=i 

Typical examples of self-affine fractals are "fern-like" and "tree-like" sets. It is rather 
difficult to find general results on the dimension of self-affine sets, and formulae such 
as (|4.6|) or ( |4.8|) are no longer valid ||100|| . From our perspective, it is sufficient to 



note that if one can show that a set J 7 is self-similar under a certain sequence of 
similarities, then eq. (|4.8| ) allows one to find its Hausdorff dimension, while if one 
only knows that a set J 7 is self-affine, the calculation of the dimension may become 
less clear. 

In the trivial example of J 77 = R D , one knows d^J 77 ) from a direct calculation 
of volumes. Then, one notices that J 7 is self-affine, but in particular it is invariant 
under certain similarities Si. Using this last property, we have recalculated the 
Hausdorff dimension. This exercise would be pointless were it not for the insight it 
can give us for the fractional Euclidean space 8%. In fact, while the symmetries of 
fractals are used to calculate or estimate their dimension, here we can reverse the 
logic and ask what are the symmetries characterizing fractional space. We know its 
Hausdorff dimension from the operational definition (scaling law of D-ball volumes). 
Therefore, because of eqs. ( |3.57| ) and (|3.2|), if £® is self-similar, then eq. ( |4.8| ) holds 
with du = dc- Sometimes, the scaling property of the measure is erroneously taken as 
the definition of self-similarity, but we have seen that scaling alone is not sufficient to 
guarantee self-similarity. Thus, we would like to go into some detail in the symmetry 
structure of fractional spaces. 

Already in one dimension, however, we see a difference with respect to the or- 
dinary Euclidean case. Take, as before, the unit interval [0,1]. If we took two 
similarities 

<S a ,i(x) := x , S a ^{x) := (1 — A) «i + a , (4-15) 



for some A and a, by eq. (|4.8|) we would obtain the expected dimension g?h = «, but 
the resulting set would be a Cantor dust. In its first iteration, it would be the union 
of the intervals [0, A 1 /"] and [a, a + (1 — A) 1 / ], with gap a — X l l a in between. In its 
second iteration, the intervals would be further split, and so on until one obtains a 
totally disconnected set. However, this is not what we expected, i.e., a continuous 
space. What went amiss is the requirement, here ignored, of endowing sets of R with 
a fractional measure. In other words, the correct procedure is to define similarities on 
the "fractional interval" [0, l] a spanned, by definition, by the geometric coordinate 
q: 

SiGz):=A<z, S 2 (g):=(l-A)g + A. (4.16) 
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These similarities would guarantee that [0, l] a = [0, A] a U[A, l] a and that d H ([0, l] a ) = 
1, consistently with eq. (|3.44 ) in one dimension. From this, one infers that dn(£ a ) = a 



and, extending to D embedding dimensions, du(S^) = Da. Calculations where 
du = D — e correspond, in the fractal picture, to regimes with "low lacunarity," i.e., 
where fractal space is almost translation invariant ||144| . 



As a mapping on x, and assuming without loss of generality that xo = [[42 
{Si[g(x)]} 1 / a oc S a> i(x) is a linear similarity, but 

{S 2 [g(x)]}^ oc S a>2 {x) = [(1 - \)x a + Ar(l + c# ^ S a , 2 (x) 

is neither linear nor a similarity. 7 Therefore, an estimate of dn(£®) via symmetry 
arguments seems unpractical in x coordinates, while it is straightforward in geometric 
coordinates. The reason why we dwelt so long on the topic of self-similarity and self- 
affinity is that it constitutes the starting point wherefrom to attack the important 
problem of the isometry group of fractional time. The integer case strongly suggests, 
in fact, that this group be given by the affine transformations 

= Stf) := A*? + a? , (4.17) 

for some D x D matrix A£ and some vector a M . We shall continue the discussion in 



[p2]| , where it will be extended to spacetimes with Lorentzian signature. 

We wish to insist upon the characterization of fractional spaces as self-similar 
sets, and rederive the result d^ = a under yet another perspective. The concept of 
self-similar measure is fundamental not only to this purpose, but also for determining 
the spectral dimension of spacetime and for generalizing the fixed-a case to multi- 
fractional scenarios. 

Let Si be N similarities defining a self-similar set T . Suppose the strong sepa- 
ration condition holds, i.e., there exists a closed set U such that Si(U) C U for all 
i = 1, . . . , N and Si(U) PI Sj^i(U) = 0. T C U is constructed taking sequences of 
similarities and the intersection of set C4 = o ■ ■ ■ o S ik (U). If \U\ = 1, then the 
diameter of the k-th iteration set is the product of similarity ratios, |Z7fc| = . . . A^. 
Let < Qi < 1 be N probabilities (or mass ratios, or weights), such that Yliidi = 1- 
One can imagine to distribute a mass on sets Uj. by dividing it repeatedly in N 
subsets of Uk, in the ratios g± : • • • : g^. This defines a self-similar measure g with 
support J 7 , such that g{Uk) — g% x ■ ■ ■ g% k and, for all sets AC J ||143|| , 



N 



g(A) = J29iQ[Sr 1 (A)}. (4.18) 



8=1 



7 By using the mean value theorem, one can show that S Uy 2 is a contraction on any compact 



interval [x a ,X(,] for x a > Xo (compare [100, example 9.8]). However, this fact does not add much to 
the discussion in the text. 
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For Cantor sets where N = 2, eq. ( 4.18|) is said to be a binomial measure (e.g. 



[p.45|| ). The case N = +00 corresponds to infinite self-similar measures, describing 
fractals with an infinite number of similarities [ |146| , |147[ . Given a real number u, we 
define the singularity (or correlation) exponent 8(u) as the real number such that 



N 



X>?A? {, ° = 1 . (4.19) 



The correlation exponent exists and is unique, since < Xi,gi < 1. As a function of 
it, 6 is decreasing and lim u _>. ±00 6{u) = =Foo. The generalized dimensions are defined 
as 

d{u) :=P^-, u^l, (4.20) 
1 — u 

and a non-singular definition, which we do not report here, is employed for u = 1. 

Self-similar measures are associated with multi-fractal sets, where the mass is not 
equally distributed among the smaller subsets of T . Fractals characterized by just 
one dimension at all scales are special cases of multi-fractals. For nested fractals, the 
probabilities are all equal to <?j = 1/N. In all fractals with equal contracting ratios 
Aj = A, the generalized dimensions all coincide with the capacity ( |11|) of the set, 
which is also the Hausdorff dimension. In fact, from eq. ( |4.19| ), 



N— = l => d(u) = -— = d c = d R . (4.21) 

Trivially, for the unit interval [0, 1] the scaling is A = 1/N — and du = 1. 

Let us apply eq. (|4.19|) to a fractional line of fixed order a. As for the ordinary 
real line, the mass is equidistributed on all subsets, and the probabilities are still 
gi = 1/N. N is arbitrary but can be fixed to N = 2, in which case g a is a binomial 
measure. In geometric coordinates, the scaling is A = 1/N = 1/2 and equals the g^, 
so 1 = 2 1_dH implies d^ — 1 and dn = a for embedding coordinates. From the point 
of view of the fractional interval spanned by the x, the scaling is 

A = 9? = (I) " = (I V , (4-22) 



N V 2 



which is smaller than in the integer case. Then, eq. ( 4.21|) yields 1 = 2 1 



consistently: the fractional charge a is the Hausdorff dimension. 
4.3 Other properties 

After 1.-3., one could mention other properties, which are more model-dependent 
and hence fail in a number of cases. For instance, one of the early definitions of 
fractals was that their dimension (defined in some of the above ways: c?b, dn, dc, 
and so on) is non- integer and greater than its topological dimension d top . There 
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are many counterexamples where the Hausdorff dimension is integer, in some cases 
smaller than or equal to the topological dimension of the set (or of the graph or 
the trail of the map it is defined by). In all cases, g?h and d top are smaller than the 
topological dimension D of the ambient space. Fractals with dn = 2 are: for d top = 1 
and D = 2, the dragon curve, the Sierpihski curve, some plane-filling curves (Moore 
curve, Peano curve), and the boundary of the Mandelbrot set; for D > 2, Brownian 
trails (almost surely, i.e., with probability 1, and g?b = g?h; however, the graph of 
Brownian motion has d-g = g?h = 3/2 almost surely ||100|| ); for <i to p = 2 and D = 2, 



the Mandelbrot set, some Julia sets, some diamond fractals, and Pythagoras tree; for 
dtop = 3 and D — 3, the Sierpihski tetrahedron. Fractals with dn = 3 and <i t0 p = 1 
in D = 3 are box-filling curves such as the Moore, Hilbert, and Lebesgue curves. 

The fractional space £® can have non-integer dimension, but its Hausdorff di- 
mension is never greater than its topological dimension, du < D. Therefore, we can 
regard fractional space as a space-filling fractal or, more suggestively, as a fractal 
associated with a diffusion process. 

Because of the empirical nature of the definition of fractals, one can conclude 
that fractional space £® can be correctly characterized as "fractal" but, since it is a 
special case of fractal with a continuous structure, it may be better to use the less 
catchy but more specific adjective "fractional." 

4.4 Fractional measures as approximations of fractals 

A conceptually independent point of view was briefly mentioned at the beginning of 
this section, where traditional fractals were found to be approximated by fractional 
measures. We can now look into greater detail at the reason why fractional calculus, 
under certain assumptions, approximates some classes of fractals. 

A rough understanding of these approximations is actually contained in eq. ( |4.15|) . 



Depending on the value of a in S a ^i we can either make a connected construction 
(a = A 1//a ) or a Cantor- type one (a = 1 — (1 — A) 1 /"). In the first case, the first iter- 
ation yields the interval [0, A 1//Qf + (1 — A) 1 /"], smaller than the desired set by a gap 
(A 1/Q + (1-A) 1 / Q ,l]. This remainder is zero only if a = 1. Therefore, we could add at 
least another similarity S3 to fill the gap, but it does not take long to convince oneself 
that the task is impossible unless one takes an infinite number of similarities as in 
fractals with infinite-type measure. This is equivalent to a continuum approximation, 
where A can be taken arbitrarily small and A 1 / Q + (l — A) 1 /" — > 1. Conversely, the first 
iteration of a Cantor-like construction gives [0, X 1/a ] U [1 - (1 - A) 1/a , 1] , with a central 
gap (A 1//q , 1 — (1 — A) 1 / ). Sending A to zero would give the same limiting set as in 
the connected construction. In a qualitative sense, we begin to recognize that real- 
order fractional measures can be regarded as approximations of self-similar fractals 
in the limit of the similarity ratio approaching zero. While deterministic self-similar 
fractals pick a countable number of ratios Aj, infinite and random fractals accept any. 
So, this limit is associated with random fractal structures. 
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This intuition is confirmed by precise arguments; we review them from the liter- 
ature, but add new comments linking independent results. Consider a function f(x) 
in D = 1 and the convolution 



Ijr(x) = v T * f :-- 



dx Vj?(x — x')f(x') 



(4.23) 



over a set J 7 C [0,1]. The kernel vjr depends on the set. For an interval, vjr is 
simply a step function, but on a fractal it can be very complicated; for example, in 
self-similar fractals vj can be determined recursively at any given order of iteration. 



In general, it is convenient to Laplace transform the convolution ( |4.23| ) 



Up) 



dxe px I T {x) =vjr(p)f(p). 



(4.24) 



If J 7 is a self-similar set, Ipip) can be expressed iteratively as an infinite intersection 
of pre-fractals, eq. ( f4.5|) . The k-th iteration has Laplace-transformed kernel Vjr(p) = 
T[n=o9 n (.P)' f° r some functions g n . Typically, for self-similar and generalized self- 
similar sets these functions are equal and with argument g n {p) = g(p^ n ), where A 
is the self-similarity ratio of J 7 . The asymptotics of g is g(z) ~ 1 + O(z) for small 
z and g(z) ~ gx + 0(z~ l ) for large z, where the constant g\ is the first probability 



weight in the self-similar measure ( 4. IS ) [78|. Then [113, 11 1] 



lim v k T (p) = vr{p) — P a F a (\np) 

k— >+oo 



where 



a 



In A 



(4.25) 
(4.26) 



and F a is a log-periodic function | |151| of period In A: 



F a (\np + n\n\) = F a (lnp) = ^ c, exp ( 27rli ^j) > ( 4 - 27 ) 

l=-oo ^ 11 ' 

for some coefficients q. The period in p is decreasing according to a geometric series. 
Combining (g~25|) with ( P~27p , 



ujr(p) = q exp [(iwj - a) lnp] 



l=—oo 



2tt/ 

hTA 



(4.28) 



Logarithmic oscillations are a curious feature of the spectral theory on fractals. In 
fact, it is known that the heat kernel trace for a Laplacian on fractals displays 



log-oscillations in the scale In a [ 152 , 153]. Oscillatory behaviour has been found 
analytically and numerically for various fractals ||154|| - ||l57| , and eq. (|4.28 ) illustrates 
a rather universal phenomenon. This is one of the most crucial points of the physical 
scenario that will emerge in ||42|| , where it shall be given adequate space. 
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Here, we are focussed only on the relation between ( [4.28 ) and real-order fractional 
integrals. The parameter a defined in eq. (|4.25|) coincides, indeed, with the fractional 
order a determining the capacity ( 4.1|) of fractional space (remember that g-y = 1/N). 
Recognizing p~ a as the Laplace transform of the fractional weight v a (x) = x a_1 /r(a) 
and comparing eq. fl4.23| ) with ( |4.25| ) and eq. (|4.22|) with ( [4.26| ), one sees that lj is 
quite similar to I a , were it not for the non-constant contribution (4.27). The main 



idea, now, is that a fractional integral of real order represents either the averaging 
of a smooth function on a deterministic fractal, or a random fractal support. The 
average of a log-periodic function F a over the period In A is defined by 



[1/2 

b a ■= (F a (]np)) := / dz F a (In p + z In A) 

J-l/2 



(4.29) 



and depends on the details of g(pX n ). Then, 



x 



a-l 



Via)' 



and 113, 114 



{hf) 



dx' Mx-x'))f(x') = b Q I a f. 



(4.30) 



(4.31) 



Taking the average in A is tantamount to dropping all the oscillatory modes in 
eq. (|4.271) except I = 0. The ui — > limit in eq. ( |4.27|) can be regarded as a large- 
Laplace- momentum limit, so that 



it 



Re(p) 



00 h r c 



(4.32) 



This is in complete agreement with [|7|, |115|| - [[L20|| , where integrals on more general 
net fractals are shown to be approximated by the left fractional integral, with the 
order a being the Hausdorff dimension of the set. The approximation (|4.32|) is valid 
for large Laplace momenta and drops all the contributions of probability weights gi 
for i > 2 out of the Laplace transform of the measure weight. These weights can 
be included to better describe the full structure of the Borel self-similar measure q 
characterizing the fractal set J 7 JE]. References |78|, |115[ |-[ P~2"0| ] make this result clear 
by a detailed Laplace analysis, but ||1 13[ |114| 1 give a sharper physical interpretation 
of what it means to take the large p limit, p is not a Fourier momentum and |p| 3> 1 
does not correspond to cutting off large scales. Yet, it is equivalent to randomize 
the fractal structure: the oscillatory structure disappears, the average of the kernel 
corresponds to the kernel itself, and the only approximation entailed in the derivation 
above is in the evaluation of the kernel ( 4.25|) (in fact, it is not obvious that Vjr*f is a 
Cauchy sequence, i.e., that the limit k — > +oo commutes with the integration). Thus, 
fractional integrals of real order represent (or, more conservatively, are intimately 
related to) random fractals. 
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Picking up again the example ( 4.15|) , eq. ( |4.25| ) states that g\ = X a . In the double 
limit A, gi — > + , the frequencies ui all vanish (the period of F a becomes infinitely 
long), a remains finite, and v^(p) ~ b a p~ a . The limit is only formal because b a (gi) 
vanishes at arbitrarily small gi, but the main point is that the oscillatory structure is 
suppressed in the limit of arbitrarily small ratio A. In this sense, fractional integrals 
of real order are continuum approximations of net fractals. This should not be 
confused with the low-lacunarity limit a — > 1. Here a is fixed and defines the 
quantity In g\j In A in the double limit A, g\ — > 0. 

It is quite the rule that fractals are nowhere-differentiable objects. In particular, 
a truly fractal spacetime is not expected to be represented by a differentiable man- 
ifold. Fractional calculus strikes quite a rich compromise between standard fractal 
geometry, where differentiability at large is given up, and a framework where only or- 
dinary differentiability is forfeited. An example from function theory is Weierstrass' 
function: it is nowhere differentiable in the ordinary sense, yet it is differentiable 



under fractional calculus ||83| . Because of the operator d n in eq. Q2.6|) , the functional 
space on which Caputo derivatives act is the space of integer-differentiable functions, 
which seems at odds with fractal geometry. If concerned by that, one could take 
another definition of fractional derivative, differing from Caputo in the absence of 



integer differentiation inside the definition, but such that d a l = |piq| . Not much 
would change in our conceptual framework, though. 



5. Spectral dimension of space 

Spectral theory is a tool to answer the famous "Can one hear the shape of a drum?" 



question ||158|| : the asymptotic spectrum of eigenvalues of Laplacian operators de- 
fined on a set provides information on the boundary of the set. This information 
is incomplete, inasmuch as "drums" with different shapes can vibrate in the same 
way, but it is nevertheless valuable. While the Hausdorff dimension depends on the 
local structure of a fractal, the spectral dimension d$ is a local probe of its topology 



[L52| , |159|| -[163| . A rigorous and fairly general definition of ds stems from the spectral 



theory on fractals ||164||— fTfiSfl. We only sketch some aspects of this theory in section 



|5.1| ; details can be found in the references. 
5.1 Harmonic structure 



Consider the topological presentation (|4.5 ) where, now, the maps «S are obtained 



from N continuous injections /j which are not necessarily the similarities of the set. 
The subsets /^(J 7 ) are weighted in two ways: by the probabilities gi of the self-similar 
measure (|4.18|) of the fractal, now constructed with the fi, and by resistance scaling 
ratios r» appearing in the definition of the Laplacian. A Laplacian /C on a self-similar 
fractal is part of the so-called harmonic structure of the set, which is characterized 
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by gi and r;. Together, these quantities are assembled into N parameters 8 



(5.1) 



We can give an intuitive coordinate meaning of the 7$. Let K(x) be a Laplacian on a 
given set defined via the injections fi(x) = rix + const. In any subcopy i, the scaling 
of the Laplacian is determined by fC[fi(x)] = r^ JC(x); summing over all the copies 
with the appropriate weight, one gets 

k{x) = Y,9irMm\ = x>? wo*)] • 



Thus, 7j 2 is the scaling of the Laplacian on the subsets. 

The harmonic structure is said to be regular if < rj < 1 for all %. When N = 2, 
this happens if, and only if, 7™ 1 = 72 2 , where n 12 are integer numbers. Given what 
one would call the standard Laplacian on a fractal, the spectral dimension rfs(J r ) is 
the unique number satisfying the relation 



N 



(5.2) 



i=i 



similar to eq. ( f4.8|) ||152|| . For simplicity, we can identify fa with the similarities 
of the embedding picture, and the resistance ratios with the contracting ratios Aj. 
Deterministic fractals have a regular harmonic structure with rj = Aj = A for all i, 
hence from (15.21) 



N 



2 



2 In AT 



ln(iV/A) d u + 1 



(5.3) 



When \ = 1/N, d$ = du = 1. The unit interval is a trivial example where the 
spectral and Hausdorff dimension coincide. 

5.2 Diffusion 

The spectral and Hausdorff dimensions of a fractal T are related to each other by the 



dimension of a Brownian motion taking place on T |[162j , |167|| - |]17T| . The anomalous 
diffusion law on fractals is characterized by the walk dimension ||162|, |163[ 




(5.4) 



8 In [152 1 and other papers, a "renormalization constant" is introduced explicitly in the definition 
of the Laplacian and of the 7^. It renders finite the discrete definition of K, in the limit of infinite 



iteration. This constant can be reabsorbed, as explained in [166] 
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Since d$ < du for a fractal, c?w > 2. The mean-square displacement of a random 
walker is a power law in diffusion time, (r 2 (cr)) ~ <j 2 / rf w_ Processes with c?w > 2 are 
of sub-diffusion, since the diffusion speed is lower than for normal diffusion (Gaussian 
process, o?w — 2). Intuitively, on a fractal sub-diffusion is due to lacunarity (i.e., the 
presence of "many holes") and/or to the very high multiplicities of certain available 
states [|172[|. Systems with <i w < 2 are caller of super- diffusion or jump processes 



p.71| , |173 , 174 , and do not correspond to fractals. 



Thus, we realize that the information obtained from the Hausdorff dimension 
was, at best, incomplete, and we cannot decide about the "fractal" nature of frac- 
tional spaces before looking at their harmonic structure or, in other words, at their 
topology, or, in yet other terms, at the way diffusion processes take place in them. 
Conversely, and contrary to popular belief in part of the physicists community, quan- 
tum gravity models with a non-integer spectral dimension do not necessarily entail 
a "fractal" spacetime. To establish whether the latter is "fractal" or not, it is im- 
portant to compare g?h with d$. To this purpose, we review the operative definition 
of spectral dimension employed in quantum gravity. The main steps are well known 
(e.g., ||175| ), but a closer and perhaps pedantic contact with the fractal-geometry 



perspective will be crucial to avoid confusion in the interpretation and construction 
of the fractional case. 

Since we want to probe the local structure of space, 9 we can imagine to place 
a test particle in it and let it diffuse in a random walk starting at point x (index /i 
omitted) and ending at point x'. For a metric space of topological dimension D and 
Riemannian metric g^ v , this process is governed by a diffusion equation for the heat 
kernel P(x, x', a), 

(d a - V 2 x )P(x, x', a) = , P(x, x\ 0) = 6{X ~ X>) , (5.5) 

where a is diffusion time (a parameter not to be confused with physical or coordinate 
time), V 2 = g~ l ^ 2 d^{g l / 2 d^) is the Laplacian defined on the ambient space, g is the 
determinant of the metric, and 5 is the U-dimensional Dirac distribution. 

The Laplacian acts on the x dependence of P, even in cases where the am- 
bient space is not translation invariant and P(x,x',o~) ^ P(x — x',o~). In fact, 
the heat kernel is the matrix element in x representation of the operator e <jV2 , 
P(x,x',o~) = (x' \e aV2 \x) . This writing means that, given a complete orthonormal 
set of eigenf unctions of the Laplacian, 

V 2 (f k (x) = -\ k <p k (x) , ^ d D Xy/g(p k (x)ip* k (x') = 5(x-x') } (5.6) 

7- J 



9 Dimension of spacetimes is always computed in Euclidean signature, i.e., after coordinate time 
has been Wick rotated. Thus, here we refer to spaces rather than spacetimes. 
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the heat kernel can be written as 



P(x, x', a) = J2 e- aXk tp k (x)vl(x') ■ (5-7) 

k 

In eqs. ( |5.6| ) and ( |5.7| ), the sum over k is replaced by an integral if the spectrum is 
continuous. 

The relative sign between a derivative and Laplacian in eq. Q5.5|) guarantees 
causal diffusion, where the particle flows away from the point x as diffusion time 
increases. The initial condition reflects the pointwise nature of the probe. Probes of 
finite shape are in principle possible, but they would not lead to a sensible definition 
of spectral dimension in the present context (in fact, we want to investigate the local 
manifold structure of a smooth space). In general, given the initial condition <f)(x, 0) 
at a = 0, the solution of the diffusion equation is 



(j)(x,0) = J d n x' y /gP(x,x',a)(j)(x , ,0) . (5.8) 

In particular, the effect of the non-local operator e°"' v2 is a shift of the auxiliary 
variable a: e°"' v2 0(x, a) = e~ a ' 9,T (j)(x, a) = 4>(x, a — a'). 

For a smooth space with D dimensions, the solution of eq. (|5.5|) can be found 
via the Fourier transform method. Consider first flat Euclidean space. The direct 
and inverse Fourier transforms of a function f(x) are 

i r+oo 

/(*) = —it / d D x f{x) =: F[f(x)} , (5.9a) 
(2ir) 2 J-oo 

i r+oo 

f( x ) = -~R d D kf(k)e lk *. (5.9b) 

(27T) 2 J_oo 

In particular, this definition is compatible with the definition of the Dirac distribu- 
tion: 

6(x) = [ d D ke ik - x , (27t)tf[(J(x)] = f d D x 5(x) e~ ihx = 1. (5.10) 



(27T) 



Phases <pk(x) = e lh ' x are eigenfunctions of the ordinary Laplacian V 2 = = 
d\ + • • • + dp with a continuum of eigenvalues, 

VV fe - = -X k e lk - X , \ k = k 2 := = kl + --- + k 2 D , (5.11) 

which allows us to write the solution of the diffusion equation as 

P{x, x', a) = J d D k P(k, a) e %Hx - x>) . (5.12) 
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In particular, due to the composition law (pk(x)(f* k (x') = e tk -( x - x ')^ the heat kernel 
( |5.7|) only depends on the relative distance of the initial and end point. P(k, a) must 
obey 



(d c + k 2 )P{k, (X) = 



lim P(k, a) 

<r-»0 



yielding 



P(k, a) = e~ ak . 
Equation ( 5.12j ) is simply a Gaussian in D dimensions, 



(5.13) 



P(x, x', a) 



(Ana 



(5.14) 



For a metric space of topological dimension D, the solution at small a is Weyl's 
expansion 



P(x, x', a) 



(Ana 



+oo 
n=l 



a 



(5.15) 



where the 0(a) remainder and the distance A(x,x') between the two points depend 
on the metric. The precise form of the coefficients A n in the expansion (sometimes 
called Hadamard-Minakshisundaram-DeWitt-Seeley coefficients) can be obtain via 
the theory of Green's functions in Riemannian manifolds [176 (for reviews, see 

The spatial average of the heat kernel P(x, x, a) at coincident points x = x' is 
the return probability 



V ^ := W) j d D x^P(x,x,a), V (D) := j d D x^, 



(5.16) 



which is the trace per unit volume of the operator e crV2 . Since the heat kernel 
P(x,x,cr) = (47ra)~ D / 2 [l + 0(a)] associated with fl5.15|) is constant in x due to 
translation invariance, to leading order in a one has V(a) ~ a~ D l 2 and the topological 
dimension of space is given by D = —2 lim CT ^ d In V(a) /din a. This formula suggests 
an operational definition, for any metric space allowing a random walk process, of 
the spectral dimension: 



(5.17) 



which can be shown to coincide with the more abstract prescription ( |5.2|) . Due to 
curvature effects, the spectral dimension changes with the diffusion parameter a, 
and one can find the phenomenon of dimensional flow at different scales already 
in classical gravity pl| . In practice, however, the solution to eq. ( |5.5| ) is difficult 
to compute for arbitrary a, and one confines her/his interest to the heat kernel 
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expansion ( 5.15 ). This expansion is valid also for spaces with boundaries or non- 
trivial topologies, in which case a must not be taken "too large" (lest global boundary 
or topology effects vitiate the estimate of d s , which must be local). On the other 
hand, random walks and Laplacians can be defined even on very non-trivial sets 
such as fractals. The exponent of the first term in Weyl's expansion is the spectral 
dimension of the fractal and it can be non- integer |]L52|| . Information on the various 
dimensions can be obtained directly from the scaling property of the heat kernel. In 
fact, under a coordinate and external-time dilation, 

P{X 2 / dw x, \ 2 / dw x', AV) = \- ds P{x, x\ a) . (5.18) 



Comparison with eq. ( |5.14 ) shows that the walk and spectral dimensions of flat space 



are, respectively, c?w = 2 and ds = D. 

To summarize, the spectral dimension ( [5.17] ) is obtained with the following in- 
gredients: 

1. An invertible transform between configuration and momentum space. In flat 
Euclidean space, the Fourier transform ( |5.9|) is a superposition of phases. 

2. A Laplacian /C. There is no unique definition of Laplacian, but the most 
natural one is such that the expansion basis of the invertible transform is made 
of eigenf unctions of /C. In flat Euclidean space, the natural Laplacian is the 
second-order operator V 2 = d^d^. 

3. A diffusion equation (T>P — K, x )P(x,x',a) = 0, which is defined by K (acting 
on x) and by the choice of diffusion process (i.e., of the derivative operator T>^, 
unrelated from eq. (|2.71|) , and of the relative coefficient between this and /C). 
These two independent ingredients correspond to the definition of a harmonic 
structure. In flat Euclidean space, diffusion is normal {V^ = d a , eq. ( |5.5| )). 

5.3 The case of fractional space 

The determination of the spectral dimension of fractional spaces highlights once again 
the main difference between fractal geometry and ordinary space constructions in field 
theory. In the first case, geometry and topology are defined by the symmetry and 
harmonic structures of the fractal, which are given at the outset. On the other hand, 
in field theory the harmonic structure stems from physical considerations. More 
specifically, we can say that the symmetry structure is first dictated by the action 
measure g and then imposed on the Lagrangian density C, but the harmonic structure 
is determined both by the symmetries and by the form of the kinetic operator. For 
instance, in ordinary field theory the natural Lorentz-invariant Beltrami-Laplace 
operator is K\ — □ = d^d^ 1 . However, any other operator of the form (/Ci) n = D n 
respects the same symmetry group, but its harmonic structure is different. Physical 
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requirements (in particular, the absence of ghosts) eventually single out K,\ as the 
kinetic operator, but this extra input is not always readily available. 

In fractional spaces, the possibility to choose a different symmetry for £ is the 
first source of ambiguity A second source is the non-unique way the diffusion equa- 
tion is defined with fractional calculus. After this preamble, let us examine the three 
ingredients of the previous section. 

5.3.1 Laplacian 



As we shall see in [|2|], symmetry arguments select two inequivalent classes of frac- 



tional field theories: 

• Fractional symmetry scenario. In ordinary integer models, the symmetries 
sym Q and sym £ of, respectively, the measure and the Lagrangian density are 
the same. If we impose sym Q = sym £ also in fractional theories, we obtain an 
action invariant under what will turn out to be the fractional generalization 
of rotation/Lorentz transformations. In this case, the invariant Laplacians we 
shall consider are 

K a := tiTVffi , K a := ^ . (5.19) 

Under a scaling transformation x —> Xx, K a — > A~ 2a /C Q and K a — > \~ 2a K, a . 

• Fractional/integer symmetry scenario. While the symmetry of the measure 
guarantees protection against proliferation of arbitrary measure operators in 
the renormalization group flow, one can prescribe a constant symmetry for C 
along the flow. Since the Lagrangian should be Lorentz invariant in the in- 
frared, we assume sym £ = sym a=1 , the integer Lorentz group in D dimensions. 
Then, 

/d :=S^d^d u . (5.20) 

Under a scaling transformation, /Ci — > \~ 2 JCi. In preparation of building the 
action functional for a scalar field in fractional space(time) [|42[| , we notice 
that the presence of a non-trivial measure weight factor 

makes the kinetic terms 0/Ci0 and — <9 M 0<9 M inequivalent. Upon integration by 
parts, the latter corresponds to 4>JC(j), where 



(5.22) 
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Several physical considerations [42, 44] will select /C, over the other choices ( |5.19| ) 
and (|5.20|) , as the Laplacian on fractional spaces. Here we do so also for technical 
reasons, which are going to become apparent to the reader: K, allows one the simplest 
and most natural analytic expansion of the fractional heat kernel in the basis of the 
fractional transform. 



5.3.2 Fractional Bessel transform 

The construction of a field theory on a fractal is subordinate to the possibility to 
define a transform thereon and to move at will from configuration to momentum 
space and vice versa. Shifting point of view on the same issue, the aim is to define a 
Laplacian operator on a fractal and study its spectrum. Indeed, the spectral theory 
of Laplacians on fractals is a hot and evolving topic in the mathematical literature. 
Currently, the state of the art is that no consistent measure-theoretic definition of a 
Fourier transform on fractals has been found, although some authors have gone so far 
as to obtain a Plancherel formula ||192||. Yet, the spectral theory has been formulated 



in a number of special cases, and general results exist for post-critically finite fractals, 
i.e., fractals which would become totally disconnected after the removal of a finite set 
of points (an example is the Sierpihski gasket; textbook introductions are ||165| , |166|| ). 



These fractals admit an actual Fourier transform [ p.93j , |194| . A natural conjecture is to 



identify the dimension of momentum space with the spectral dimension ds, calculated 
by looking at the spectral properties of the Laplacian. This can be motivated, for 
instance, in thermodynamical systems, where the spectral dimension comes from the 
momentum space trace needed to implement thermodynamical equations of state 
195| |. In general, ds differs from the Hausdorff dimension g?h of configuration space. 
On a more phenomenological ground, there are some results also in particular 



models of spaces with fractal-like features. In ||196|| , an invertible Fourier transform 
was defined on a metric space with non-integer dimension, only for the class of 
function generated by the Gaussian and used in perturbative field theory. On the 
same class of functions, an invertible transform exists on spaces with translation- 



invariant Lebesgue-Stieltjes measure [Q. Then, the measure in momentum space 



has the same dimension as in configuration space. This agrees with [|39|| , where we 
defined a transform with Stieltjes measure such that the engineering dimension of the 
measure is the same in both momentum and configuration space and ds = g?h- We do 
not corroborate the same conclusion in this paper, in fractional spaces, transport can 
be anomalous with respect to the natural differential structure and, by the diffusion 
equation method, we shall find that ds = du only in special cases, but not generally. 
In these spaces, the spectral dimension is not the dimension of momentum space. 

The derivation of an invertible transform in fractional spaces considerably differs 
from that of [|39] in the final result. The discrepancy lies in the fact that a naive 



replacement of the measure in the Fourier integral does not lead to an invertible 
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transform, 10 due to the loss of translation invariance in the measure (such an invari- 
ance was formally assumed in ||99|| ). We give here the correct answer, which is proven 
in a separate publication pi |. One can show that, given the fractional measure 



dg a (x) = d xv a {x) 



(5.23) 



with weight (|5.21|) , T a (k) = g a (k) and the fractional transform and anti-transform 
are 



/>) := 

m - 



dg a (x) f(x)c a (k,x) =: F a [f(x)} 
dg a (k) f(k) c a (k,x) , 



where 



- a (k,x) := \{c a ^kx) := J]r(a)(A;V) 1 -f J±- X (k»x») 



(5.24a) 
(5.24b) 

(5.25) 



and J v is the Bessel function of the first kind of order v. Notice that, consistently, 
Ci(k,x) = (2/n) D / 2 cos(A;ix 1 ) ■ • ■ cos(k D xo), and one recovers the ordinary Fourier 
cosine transform. 

The check that the fractional Bessel transform thus defined is invertible makes 
use of the integral representation of the ordinary Dirac distribution in terms of Bessel 
functions [44]. In one dimension, 

r+oo 

5(x — x') = x / dk k.J v {kx)J v {kx') , (5.26) 
Jo 

for Re(^) > — 1 (a > in the above calculation). In particular, in D dimensions and 
for v = a/2 — 1, eq. (|5.26| ) yield the correct integral representation of the "fractional" 
Dirac distribution: 



5 a (x, x) := v a 1 (x)5(x — x) = / dg a (k) c a {k,x)c a {k,x') , (5-27) 
which is not translation invariant. From eqs. ( |5.24b| ) and ( |5.26| ) , in one dimension 

yo 

dg a (x) S a (x, x')[x a/2 f(x)} = x' a/2 f(x') , (5.28) 



and the generalization to D dimensions is straightforward {x a l 2 — > (x 1 ■ ■ ■ x D ) a ^ 2 , 
and so on). 



10 



I am indebted to G. Nardclli for pointing it out. 
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5.3.3 Diffusion equation 

Depending on the choice of Laplacian and of the operator the diffusion equation 
will describe different diffusion processes, characterized by different values of the 
walk dimension. When anomalous dimensions are involved, there is no unique way 
to determine the diffusion equation, except on heuristic or phenomenological grounds. 
Anomalous diffusion can be realized also without fractional derivatives (e.g., [ 125 |), 
but here we pick the form 

{pl -K,)P{x,x',a) = 0, P(x,x',0) = 5 a {x,x'), < < 1 , (5.29) 

where 

P^R-toof}, (5.30) 

and we chose the fractional delta as the natural initial condition. Generalized to 
{T>P — IC.y)P(x,x' ,a) = 0, eq. (|5.29|) is often called fractional wave equation. When 
(3 = 1 = 7, this is ordinary diffusion (Brownian motion). When /3 ^ 1 and 7 = 1, 
it corresponds to fractional Brownian motion. Finally, when = 1 and 7 ^ 1 the 



process is a Levy flight [134 



5.3.4 Spectral and walk dimensions 

We have all the ingredients to determine the spectral dimension of fractional space 
with the heat-kernel method. The geometry of a set affects the spectral dimension 
in two different ways: through the fractal structure and the metric structure. In a 
fractional theory, what we call fractal structure is actually the calculus associated 
with the fractional manifold, and the metric structure is given by the fractional 
metric. Drawing inspiration from the ordinary manifold result, we can ignore curva- 
ture effects to obtain the spectral dimension of fractional spaces for small diffusion 
parameter. 

Solutions of the diffusion equation ( |5.29|) and its generalization (T>f — JC~ f )P = 



are known for phenomenological models of statistical mechanics, where an anomalous 



diffusion process takes place in ordinary space |46|, [71, [131|| - fi"3"4[| . Here, the very 
geometry of space is modified, and we need to start from scratch. 

Given a diffusion equation (integer or fractional), a most practical way to solve it 
is to render it algebraic. In ordinary space, one writes the heat kernel as a quadratic 
superposition of eigenstates (pk(x) of the Laplacian, eqs. ( |5.6| ) and (|5.7|). By con- 
struction, the Laplacian is chosen so that its eigenfunctions (phases, cosines or sines) 
also constitute the basis of the momentum expansion of the Fourier transform. 

In the fractional case, the eigenfunctions of the Laplacian K\ are not the functions 
c Q (k,x) in the Bessel transform ( 5.24a|) , so one is left with two different options: to 



expand the heat kernel in the Laplacian eigenstates (fk{x) or in the functions c a (k, x). 
The first possibility seems the one guaranteeing an algebraic solution of the diffusion 
equation, but in fact it is not compatible with the initial condition in ( |5.29| ). We 
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omit the proof that the expansion (|5.7| ) does not work properly. The second case 
turns out to be both analytically exact and compatible with the initial condition, but 
only for the Laplacian /C. The c a are eigenfunctions of the kinetic operator ( |5.22| ), 



)Cc a (k, x) = —k 2 c a (k,x) . (5.31) 

This is a Bessel equation with two independent solutions given, in one dimension, by 
(kx)~ u J v and (kx)~ u Y lJ , where v = a/2 — 1 and Y v is the Bessel function of the second 
kind, or by their complex linear combination into Hankel functions. Transforms 
can be defined also via Y and the Hankel functions, but only the Bessel transform 
with J is invertible and unitary, i.e., it is the only one where the Dirac distribution 
admits a Bessel integral representation. Therefore, the existence of an invertible 
unitary transformation between configuration and momentum space determines the 
natural Laplacian in configuration space. Notice that the c a already have the correct 
normalization (eq. ( |5.27| )). 

The heat kernel ( |5.7| ) is generalized to fractional space by identifying <fk( x ) with 
c a (k,x), by replacing the sum over k with an integral with the momentum-space 
measure g a (k), and by replacing the exponential with the eigenfunction fk(a) of the 
operator (|5.30|) , with eigenvalue —k 2 : 

(V^ + k 2 )f k (a) = 0. (5.32) 

For different choices of T>%, the solution is 



h{°) = { e-*-*) 1 " , V? = J% , (5.33) 

all in agreement with eq. (|5.13|) when /3 — y 1. The final form of the fractional heat 
kernel is then 

r+oo 

P(x,x',a)= dg a (k) f k (a) c a (k, x)c a (k, x') . (5.34) 
Jo 



By construction, it obeys the diffusion equation ( |5.29|) and is not translation invari- 



ant. The initial condition is also correct, as one can see by comparing P(x,x',a) 
with eq. (|5.27 ). Furthermore, the general solution $(x,a) of the diffusion equation 



for the initial condition $(x, 0) is always of the form 

*(x,a)=(j[ /i x' i y<f>(x,a) ) (5.35) 
and is naturally given by the fractional generalization of eq. (|5.8|) , 

r+oo 

$(x,a)= dQ a (x')P(x,x\a)§(x'^), (5.36) 
Jo 
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in accordance with eq. Q5.28 ). 

The return probability is the fractional counterpart of eq. (|5.16|) and reads 



n 



VaJO 



dg a (x) P(x, x, a) 



+00 



v , dg a (x) 

Jo JO 



dQ a (k) f k (a) c 2 a (k,x 



(5.37) 



where V a := J dg a (x) is a divergent total volume prefactor. For each of the D 
directions, one has to solve a double integral of the form 



dxx dk k fk(a)J„(kx) . 
Jo 

In the ordinary flat case, the coordinate dependence factorizes and is cancelled out by 
the volume prefactor, but here the integrals in x and k do not commute. The choice 
of integration ordering and the regularization of eventual divergent contributions do 
not affect the spectral dimension, eq. ( |5.17| ), because from our choices of fk{&) one 
can always rewrite eq. ( |5.37| ) as a function of a times a pure (a-dependent) number. 
Let (3 include also the value 1. From eq. ( |5.33|) , upon the rescaling 



k 



a^k. 



x 



a 2 x . 



one can always factorize eq. (|5.37|) as 



V(a) 



DaB 

a 2 



1 



d D x / d D kC{~k,x;a) 
Jo 



(5.38) 



(5.39) 



The a dependence comes exclusively from the volume prefactor, while the numerator 
of (|5.37[) is rendered dimensionless by the rescaling ( |5.38|) . Integration of the a- 
dependent function C entails an ordering and, eventually, a regularization choice. 

Combining eq. ( 5.39 ) with ( |5.17 ), we finally obtain the spectral dimension: d$ = 
Da(3 for = d a , and ds = Dotj3 for T>P = 008%, d%. In general form, 



ds = I3d u . 



(5.40) 



Some remarks are in order: 

• For normal diffusion (/3 = 1), 



(is = di 



Da 



(5.41) 



realized by integer operators d a and K. The order of the diffusion operator is 
half that of the Laplacian. 
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• When /3 < 1, diffusion is anomalous but fractional space can be regarded as a 
fractal. 

• The case f3 > 1 does not correspond to a fractal, since ds > d^, the operator T>^ 
is higher order from the point of view of the differential structure of fractional 
space, and it is responsible for super-diffusion. In particular, for /3 = I /a the 
spectral dimension coincides with the topological dimension of space. 

In all these cases, the spectral dimension is constant and non-vanishing. There 



is no contradiction with the findings of ||197 |, where the spectral dimension does 
not converge to the embedding dimension when the lacunarity of fractals becomes 
asymptotically zero. Here, fractional spaces are not low-lacunarity approximations 
of fractals and we do not take the limit A — > 0. 

The results are summarized in table ^ for ds and the walk dimension g?w 



VP 


ds 


dw 




dn 

(3d K 


2 

2//3 



Table 4: Spectral dimension ds and walk dimension dw of fractional space £® for the 
natural Laplacian ( 5.22| ) and different diffusion equations. Fractional space is fractal only 
if d w > 2 03 < 1. 



6. Discussion 

The results presented in this paper are the first step towards a field theory on frac- 
tional spacetimes. We have detailed the geometric properties of the fractional equiva- 
lent of Euclidean space, with no time, no matter, no gravity, and fixed real dimension. 
The goal, of course, will be to include the physics in a controlled way. Bits of it, such 
as the generalization to spacetimes with Lorentzian signature, only require minor 
modifications of the fractional construction. Others, such as the formulation of a 
multi-fractal scenario, the recovery of four dimensions at large scales, and inclusion 
of the log-oscillations of fractal geometry, are not difficult but entail a major change 



of perspective. We endeavor to complete this programme in |[42|| , still in the absence 
of gravity. 

We would like to conclude with a remark at the interface between quantum 
gravity and mathematics. It has been recognized that effective spacetime emerging 
from quantum gravity scenarios has a scale-dependent spectral dimension, a feature 
on which part of the community has grown the belief that quantum spacetime is, 
somehow, "fractal." However, aside from the spectral dimension, fractal geometry 
is an arsenal of tools which has been scantly exploited in quantum gravity. As a 
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consequence, a deeper understanding of fractal properties of spacetime has seldom 
gone beyond qualitative remarks based on quantitative determinations of ds- Taking 
advantage of this arsenal in a less frugal way would open up a wealth of possibilities, 
as we shall argue in the companion paper. On the other hand, progress in pure fractal 
geometry is very much ongoing and the effective insights in the physics literature can 
suggests mathematicians some interesting directions of research. For instance, to the 
best of our knowledge there is no systematic formulation of a Fourier transform on 
fractals. To a physicist, transforming to momentum space is important both for 
doing field theory and for computing the spectral dimension. Related to that, we are 
unaware of any good physical transport model where jump processes naturally occur. 
These are non-local diffusing processes characterized by discrete jumps, rather than 
continuous movements. While for local diffusion the walk dimension is bounded by 
2 < dw < d}i + 1, in non-local diffusion < c/w < <^h + 1 ||171|| . The cases where 



fractional spaces have <i\y = 2g?h/g?s < 2 might correspond to continuum models of 
jump processes, but the physical meaning of this is presently unclear. Yet, transient 
regimes where d s > <i H do arise in other approaches to quantum gravity, as in causal 
dynamical triangulations |16| or in non-commutative spaces |[L98|| . Further study of 



the subject promises to be stimulating. 
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